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We give an extensive treatment of the Constant Mean Curvature (CMC) Ein- 
stein flow from the point of view of the Bel-Robinson energies. The article, 
in particular, stresses on estimates showing how the Bel-Robinson energies 
and the volume of the evolving states control intrinsically the flow along evo- 
lution. The treatment is for flows over compact three-manifolds of arbitrary 
topological type, although the form of the estimates may vary depending on 
the Yamabe invariant of the manifold. We end up showing well posedness of 
the CMC Einstein flow with H'' x regularity, and proving a criteria for a 
flow to be a long-time flow on manifolds with non-positive Yamabe invariant 
in terms only of the first order Bel-Robinson energy. 
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1 Introduction. 



The aim of this article is to provide an intrinsic treatment of the Constant Mean 
Curvature (or simply CMC) gauge, entirely in terms of the space-time Bel-Robinson 
energy and the space-like volume. Roughly speaking, the CMC gauge foliates the 
space-time in such a way that, at every leaf, the local rate of (normal) volume ex- 
pansion is constant (i.e. independent of the point). The space- like volume, in turn, 
is linked to the space-like scalar curvature and (though it) to the Yamabe invariant 
of the three- manifold. As it turns out, the link between the space- like volume and 
the spatial geometry is strengthened if the Bel-Robinson energy associated with 
the Riemann tensor of the space-time is taken into the picture. We will display 
precise estimates showing how the volume and the Bel-Robinson energies control 
the geometry of the states. This, in particular, makes the Bel-Robinson energies 
and the space-like volume a set of appealing variables to control the CMC flow 
along evolution. We will give a detailed discussion of intrinsic as well as elliptic 
estimates. 

The CMC gauge and Bel-Robinson energies have been used together several 
times in the past. In [13], for instance, Christodoulou and Klainerman approached 
the stability of the Minkowski space-time using a maximal foliation and a elaborated 
control of the Bel-Robinson energies of appropriate Weyl fields. In their work, Weyl 
fields were used in particular as the fundamental variables from which to reconstruct 
the space-time metric. On the other hand, in [4J, Andersson and Moncrief gave a 
proof of the stability of flat cones following essentially the same argumental lines. 
This case, which can be considered as a compact version of [13], is however greatly 
simplified thanks to the expansion in volume (and the compactness of the CMC 
Cauchy slices). In the context of the initial value formulation of the Einstein 
theory, Weyl fields were used by Friedrich in [T^] • Friedrich included the conformal 
Weyl tensor of the space-time as a variable, and by doing so he obtained different 
hyperbolic reductions of the Einstein equations from which to launch initial value 
formulations. We will take, in spirit, several elements from these works. Namely, 
we will study how to control the space-time metric from the Bel-Robinson energies 
(of suitable Weyl fields) and use that knowledge to give a treatment of the Cauchy 
problem in General Relativity entirely inside the framework of Weyl fields. 

To step further in the description of the contents, let us introduce some ter- 
minology first and then state, to exemplify, some of the main estimates that will 
proved. The reader can consult the background section for a detailed account on 
notions such as Bel-Robinson energy, or harmonic radius that we will mention be- 
low. Consider a cosmological solution (M,g) and a Cauchy slice E. Denote by g 
the spatial three-metric and by K the second fundamental form of S. As is well 
known K — —\g' where g' is the time derivative of g in the normal direction to 
E. Thus it will be natural to call the pair {g,K) — (3,-^5) a state. A state is 
in particular a CMC state if the constant mean curvature k — trgK is constant on 
E. Consider the space-time Riemann tensor Rm. If the space-time solution g is a 
vacuum solution the curvature Rm is a Weyl field as are its covariant derivatives 
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V^Rm in the normal direction T to E. Consider their associated Bel-Robinson 
energies and denote them by Qj . One of the main results will be to show 

Lemma 1 (Sobolev norms vs. Bel-Robinson norms) Say I > and say S is 
a compact three-manifold. Then the functional (defined over states {g,K) with 
k = -3; 

1 1 

||(g,X)||sR = - +V + ^Q, = - + V + Ej, 

^ j=o ^ 

controls the H^'^^ -harmonic radius rj^2 ^'^'^ ^^^^ H^'^^ -norm of K . 

Above, V = {^YVolg{Ti) is the so called reduced volume (see later), which is 
essentially the usual volume if k was fixed and equal to —3. v is the so called 
volume radius which is a measure of the local "flatness" (see later). H^'^^ is the 
I -\- 1-Sobolev space and the subindex g indicates that the Sobolev-norm is the 
natural constructed out of the metric g (see the background section). The proof 
of this Lemma emerges basically as a corollary of a series of partial results which 
have, however, independent value. Consider for instance the case / = in Lemma 
[TJ We show in Proposition [6] that the iJ^-norm of K, namely, 

\\Krm= I \^K\' + \K\'dvg, 

is controlled only by V and Qo- This will follow from the explicit intrinsic estimate 

/ 2|Vif|2 + |X|4dwg<C(|fc|V-l-Qo), 

in the case the Yamabe invariant Y{Yj) is Y{Yj) > 0, and 

/ 2\vk\'' + \K\^dvg < C{\k\iV - V„/) + Qo), 

if Y{Yj) < 0, where K is the traceless part of K, C is a numeric constant and Vi„/ 
is the infimum of the reduced volume V among all CMC states {g, K). It turns out 
that Vinf = (^^^(S))2 (see the background section). Being formal, an intrinsic 
estimate is an inequality (or equality) involving intrinsic Sobolev norms H* but not 
requiring the volume radius v. The intrinsic estimates on K above, imply, as we 
prove in Proposition [51 an intrinsic estimate on the L^-norm of the Ricci curvature 
of the three- metric g. Namely we show 

\\Ric\\l2 <Ci\k\V + Qo). 

It follows from the fundamental theorem of convergence of Riemannian manifolds 
(see the background section) that a priori control on V and Qq, and in addition a 
priori control on v, implies control on the iJ'^-harmonic radius r2. This is the way 
Lemma[I]is proved when / = 0. The situation when J > requires the use of elliptic 
estimates and therefore a priori control on v. For instance when I — 1, it is true that 
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V, V and E\ control ||fir|l//2 but through an inequahty that involves implicitly the 
use of and which arises through the use of suitable elliptic estimates. To handle 
elliptic estimates when the background metric g is one of the variables, we use 
Sobolev norms defined with respect to atlas composed of harmonic coordinates. 
Being more precise, we define W- canonic harmonic atlas, and measure Sobolev 
norms with respect to them. A chart {xa } in a iJ*-canonic harmonic atlas, has the 
property that if we scale g as g = and {xa\ as = then Xa gets defined 
over a ball of radius one in the metric g and moreover the usual iJ|^-norm of g is 
bounded above by a fixed (but arbitrary) constant. Supplied with some additional 
technical requirements, _ff*-canonic harmonic atlas provide us with a more or less 
standardized way to define Sobolev norms with respect to atlas made of harmonic 
coordinates. We will denote by Hj^ Sobolev spaces defined with respect to an 
atlas A and in particular with respect to a canonic harmonic atlas A. Once an 
estimate, or a certain inequality is proved using the Sobolev spaces H^, one can 
guarantee that a similar inequality or estimate holds with respect to the intrinsic 
norms iJ*, but the given estimate or inequality involves necessarily v. For instance, 
in Proposition [14] we prove that V and Qq control the 77|-norm of 1/N where N 
is the lapse function in the CMC gauge when one take fc as a choice of time. This 
(not self-evident) result, which in particular implies that N is never zero (even for 
states (g, K) with x regularity) and controlled from below by V and Qo, 
requires the intermediate use of the norms i/J^, where ^ is a _ff ^-canonic harmonic 
atlas. 

Lemma [1] naturally points towards a proof of well-posedness of the CMC Ein- 
stein flow, entirely in terms of the intrinsic quantities V and £j. In this respect 
we will prove the following version of well-posedness of the Einstein CMC-flow. 

Theorem 1 Let (S,-4) he a H'^ -three-manifold and say (gQ^Ko) is an initial state 
in X H"^ with <0. Then 

1. There is a unique H^-flow solution over an interval I = (fc-i, fci) with — oo < 
k-i < ko < ki < 0. Moreover the size inf{\k-i — fco|, |A;o — fci|} of the time 
interval on which the solution is guaranteed to exist is controlled from below 
by l/v{ko), ln|/co|, V{go,Ko) and Siiko). 

2. There is continuity with respect to the initial conditions if we measure the 
space of initial conditions with the x norm and the space of flow 
solutions with the BR-norm. 

3. Because of item 1 above, we have the following continuity principle: a flow 
solution {(g,K), (N, X))(k) is defined until past of k* < (or before fc* < «/ 
the flow is running in the past direction) iff limsupk^k^l/h' + V{k) -\-£i{k) < 

CO. 

The BR-norm (Bel-Robinson norm) of a flow in item 2 of the Theorem above is 
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defined by 

fe=i 

(1) l|g||s_R = ll.9llco(/,2)(ff^) + \\K\\co{IS){Ha) Bk)\\co{I,0){HA)^ 

where {Eq, Bq) and (£'1, Bi) are the electric and magnetic components of the Weyl 
tensors Wq — Rm and Wi respectively. The atlas A used in ^ is the atlas that 
is given from the manifold {Tj,A). The space C{I,j){Hj\) is defined as usual as 
the space of continuous functions over an interval / with values in the i/^-Sobolev 
space of tensors of a fixed rank. 

The initial value formulation of General Relativity in the CMC gauge has been 
considered in the literature some times in the past. In particular in |12j3, Choquet- 
Bruhat and York gave an exposition of the Cauchy problem applicable to the CMC 
gauge and for initial states with x regularity. In [5], Andersson and Mon- 
crief considered the CMC gauge over hyperbolic-three manifolds, in the spatially 
harmonic gauge and for initial states with i/2.5+e ^ regularity. Theorem [1] 

partially (but not entirely) overlaps the results in [T^] and [S] . 

Note from ^ that the flow of three-metrics g{k), is measured in and not 
in H^. In our approach we cannot guarantee that the metrics g{k) will lie in 
(see however the claims in |12j). To circumvent this problem we included the terms 
{Eo,Bo) and {Ei,Bi) in Indeed, by LemmalU ^ + V + Qo + Qi controls the 
r3-harmonic radius and therefore (see the background section) the -norm of 
g where Ah is a iJ "^-canonic harmonic atlas. As ||5||//2^ controls ^ and V (when 
|fc| is bounded) it follows that the "norm" ||g|| ^^2^ -I- || (£^0, 5o)||l^ + II (^^i, ^i)||l^ 
controls the norm || (7II//3 .To transform this fact into a technical tool we introduce 

in Section [33] a dynamical iJ "^-harmonic atlas A{k). The BR- norm acquires then 
special relevance when A{k) is taken into account. The dynamical atlas will be a 
fundamental piece in all our treatment of the Einstein CMC flow. 

Consider a CMC initial state {go,KQ) with ko ^ over a three-manifold 
that we label as Sq. It is well known the state {g, k) gives rise to a C°° space-time 
(M, g). By the method of barriers it is showrH in the background section that there 
is a unique C°° CMC foliation of M around the slice Eq where k varies strictly 
monotonically. Thus, for C°° space times arising from initial C°° CMC states, the 
CMC gauge is well defined at least on a neighborhood of the initial slice. The 
strategy we will follow to prove Theorem [1] is to show that the completion of the 
space of C°° CMC flow solutions under the BR-norm ([T]) is precisely the space 
of H^-CMC flow solutions (see Deflnition [4] in Section [3^ . Note that this is not 
the traditional approach to prove well-posedness of a PDE but it illustrates the 
applicability of the estimates in Section [31 

There is an underlying reason of why to describe the evolution in terms of the 
intrinsic quantities V, v and £1 and this has to do with the continuity principle as 
in item 3 of Theorem [T] A common strategy to analyze the long-time evolution of 

■^We would like to thank the referee for pointing out this reference to us. 
•^This fact is well known (see the discussion in [23)). 
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flows is to design a suitable continuity principle in terms of appropriate "continuity" 
variables that could "describe"]^ the flow when at least one of the variables breaks 
down. In other words, and being slightly ambiguous, if one defines the notion of 
singularity as equivalent to the blow up of the "continuity" variables, then no extra 
information other than the one provided by the variables themselves, would be 
required to analyze the structure of the singularities. In still vague but heuristic 
terms, one would like to have a "complete set of continuity variables" . The right 
choice of the "continuity" variables is thus of central importance. Indeed this is 
one of the major obstacles to understand the long-time evolution of the Einstein 
flow. Although it is unlikely that the quantities v, V and Ej are best adapted for 
this purpose, we believe they conform a interesting set, deserving to be considered 
with care and in depth. A related problem of considerably technical as well as 
conceptual difficulty, is whether £q = Qq can be used instead of £i = Qq + Qi 
in Theorem [TJ This problem is known (with variations in the formulation) as the 
-ff^-conjecture. Its difficulty lies however well beyond the scope and the techniques 
developed in this article. Still some of the estimates here developed may be of 
tangential interest in this goal. 

We require, as part of the definition of flow solution (Definition |4] in Section 
I3.4|l . that a {H^} flow {g,N,X) (of three-metrics, lapse and shift) is in fact the 
flow induced by a CMC foliation on a {H^) space-time (M,g) (see Definition [3] in 
Section [5^ . It follows from Theorem [T] that an initial data {gQ,Ko) in x 
gives rise to a unique flow solution in a given shift X or, equivalently, a unique 
space-time (M, g) up to space-time diffeomorphisms. We note that the equivalence 
between flow and space-time solutions is not always imposed as a requirement in 
some treatments of the initial value formulation in General Relativity. 

As the variables v, V and Qo are intrinsic to the space-time, we can see that 
Theorem [T] is in fact independent of which shift X we take. Indeed, our approach 
to the Einstein flow is independent of the shift. For this reason we have introduced 
a general notion of admissible shift, made up (in some sense) out of the minimum 
requirements that a function X(g,K) must satisfy to be the shift of an Einstein 
CMC flow. We may well take X = (which is an admissible shift) all though the 
article. 

Let us give now an overview of the contents of the article. In the background 
section we provide a detailed account on the notions of: Einstein flows (Section 
12. ip . the CMC gauge (Section [221), Weyl fields and Bel- Robinson energies (Section 
12. 5p and scahng (Section [2.6p . Although most of the article is based on standard 
functional analysis, we introduce some new terminology (as canonic atlas) that 
needs to be presented with care. This is done in Section 12.41 Section [3] contains 
the core of the article. In Sections 13.11 and 13.21 we discuss the notion of reduced 
volume and give a cosmological interpretation of its monotonicity. The fact that 
the reduced volume is monotonically decreasing in the expanding direction has 

*Note that "describe" in GR means "describe the space-time" an not only the flow. This fact 
makes GR a particularly difHcult theory to treat only from the point of view of dynamical flows. 
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important consequences for flows over three-manifolds E with non-positive Yamabe 
invariant. One can read an easy implication of the monotonicity in item 3 of 
Theorem [T] In fact, because V is decreasing, it remains bounded to the future, 
and therefore we can dispense with it in the continuity criteria. In Section 13.31 
we show how the variables v, V and £i control CMC states (5, K). To clarify the 
difference between intrinsic estimates and estimates of elliptic type, we have divided 
the section into two subsections. Sections 13.3.11 and Section r3. 3. 21 In Section r3. 3. II 
we introduce the intrinsic estimates. These estimates are in terms of V and Qq 
only. In Section 13.3.21 we introduce the estimates of elliptic type. This necessitates 
of a laborious treatment of the currents J(Wi) associated to the Weyl fields W^. 
It ends up with a proof of Lemma [T] Section 13.41 is essentially an extrapolation of 
the analysis done in Section [3731 for single states (g, K) but now for flows (g, K){k). 
We introduce a BR-norm and show how it controls the flow and the space-time. In 
Section [3.41 too, we introduce the notion of a dynamical harmonic atlas. We spend 
much of Section [33] analyzing this concept. The dynamical harmonic atlas is the 
necessary ingredient to show that flow solutions are actually space-time solutions. 
Section [47T] is devoted to the initial value formulation. All the treatment is entirely 
in terms of Bel- Robinson energies and Weyl fields. Finally, we present in Section 
14.21 some partial results on long time CMC flows over manifolds with non-positive 
Yamabe invariant. A long time CMC flow is one for which the range of k contains 
an interval of the form (/cq, ojl. What we prove is the following 

Theorem 2 Say < and (M ~ E x M, g) a smooth C°° maximally globally 

hyperbolic space-time. Say {g,K){k) with k G [a,b) is a CMC flow where b is the 
lim sup of the range of k and say £1 < A. Then if (M, g) is future geodesically 
complete the CMC flow is a long time flow i.e. 6 = 0. 

Except in a few exceptions, all the article deals with the vacuum Einstein theory. 
In those cases where matter is considered, we will do so assuming that the energy- 
momentum tensor satisfies the dominant energy condition. 

2 Background. 

2.1 The (vacuum) Einstein flow. 

All manifolds and tensors in this section and the next (Section 12. 2p are con- 
sidered to be smooth i.e. C°°. We will denote space-times (i.e. formally a four 
manifold M and a Lorentz ( — h -I-H-) metric g) by (M, g). All through, space-time 
tensors will be boldfaced. We will think the Einstein equations in vacuum 

Ric - ^Rg 0, ( or Ric = 0), 

^It is conjectured that any flow over a manifold S with non-positive Yamabe invariant is in 
effect a long time flow 1231 (see also [6] and refferences therein). 
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not as global equations on (a given) M but rather as evolution equations on mani- 
folds of the form M ~ E x E, where E is an orientable and compact three-manifolc|§. 

Say M is 3+l-split by the diffeomorphism : E x R — > M, in a way that 
'P* {s)\sx{t} = g{t) is Riemannian for every f(i:M = ExM^Mis the time 
function coming from projecting into the second factor). Writing dt = NT + X 
with T a normal to the time foliation and X tangential to it, the space-time metric 
splits into E and dt components as 

(2) 0*g = -{N"^ - \X\^)dt^ +dt®X* +X* (ddt + g. 

In the formula above g was extended to M to be zero along dt (i.e. zero when one of 
the entrances is dt) and X* — X^gab- N is called the lapse and X the shift vector. 
As seen in a 3-l-l-splitting, a space-time metric g is characterized uniquely by a, flow 
(a path) of three-metrics g{t) and lapse-shift {N,X){t) on a fixed manifold E. We 
will denote g4,{t) = {{g), {N,X)){t) and call it the Einstein flow in the (l)-splittin^ 
It is convenient to include the second fundamental form of the time foliation as part 
of the definition of the Einstein flow, thus = {{g,K), {N, X)){t). The Einstein 
flow equations in vacuum are 

(3) R=\Kf~k'^, 

(4) VK~dk^ 0, 

(5) g^-2NK + Cx9, 

(6) k = -VVN + N{Ric+kK ~2KK) + CxK. 

Equations ^ and Q are the constraint equations and equations (|3.ip and ^ the 
Hamilton- J acohi equations of the flow. Contracting ([6]) and using the constraints 
we get the lapse equation 

(7) - l\N +\K\'^N = k- X{k) = N'^rk. 

This equation can be seen too as the second variation of volume at the leaves of 
the time foliation. 



° Solutions in this kind of topology are called cosmological solutions. As usual, hypersurfaces 
difFeomorphic to S having Riemannian inherited metrics g are called Cauchy hypersurfaces 

^Prom a Lagrangian point of view one studies the flow g(t) = ({g, v), {N, X)){t) of position g, 
velocity v = g = Cg^g = —2NK + Cx9 {K is the second fundamental form of the time foliation, 
note the sign convention) and lapse-shift (A'^, X) as a solution to the Euler-Lagrange equations of 
the Einstein-Hilbert action (see 1241 ) S = /j^^j^Rdiig -I- 2 Jq^y.^.!) ^<^'"g- ^^gK = A; is the mean 
(extrinsic) curvature. Alternatively, from a Hamiltonian point of view one studies the evolution 
of the flow = {{g, P), {N, X)), of position g, momentum P = {K — kg)dvg and lapse-shift 
(A', X) as a solution to the Hamilton-Jacobi equations of the Hamiltonian (make P = {K — kg)) 

?i = Jjn A''(|Pp ~ ^ ^) ^ 2X .(\/ .P)dvg , generated by the energy and momentum constraints 

functions S = \K\^ — — R and V = —2SJ.{K — kg), which are well know to be conserved along 
the flow and identically zero (in vacuum). 



8 



It is well known that the constraint equations comprise the necessary and suf- 
ficient conditions on initial states {g, K) to have solutions to the Cauchy problem 
for the Einstein flow equations. Given a solution (M, g) of the Einstein equations 
and a diffeomorphism 0, then the flow is a solution of the Einstein flow equa- 
tions, in particular two different diffeomorphisms 4>i and (f)2 that coincide and have 
same differential over an initial Cauchy hypersurface and therefore inducing the 
same initial states {{gojKo), {No,Xq)) give rise to two different solutions g^^ and 
g02 on S X M. Their associated space-time metrics from equation ([2]) are however 
isometric. This is the only freedom, i.e. module diffeomorphisms, solutions are 
unique. Note that we can always adjust the diffeomorphism over an initial Cauchy 
hypersurface to make N and X have any prescribed values (with > \X\'^). We 
now state the well known theorem of well posedness of the C°° Cauchy problem 
(see for instance [l8],[2l] and references therein). 

Theorem 3 (Well posedness of the (C°°) Cauchy problem) (Existence) Given and 
initial state (go, Kq) satisfying the energy and momentum constraints and arbitrary 
initial lapse-shift {Nq,Xq), there is a solution to the Einstein flow equations 
(0j over a short period of a parametric time. (Uniqueness) Two globally hyperbolic 
solutions gi, g2 to the Einstein flow equations with the same initial state (goi^o) 
and same initial lapse and shift (Nq,Xq) are isometric i.e. there exists (M,g) and 
diffeomoprhisms (pi and (f)2 such that g = 0*(gi) = 0*(g2)- 

The freedom in the choice of 0, i.e. the invariance of solutions under diffeomor- 
phisms is called gauge freedom. A choice of from g in such a way that if gi is 
isometric to g2 then 0(gi) = 0(g2) is said a choice of gauge. 

2.2 The CMC gauge. 

A state (g, K) satisfying the energy and momentum constraints with k = trgK 
constant is called a CMC state. A cosmological solution on E x R is in the CMC 
gauge if the extrinsic mean curvature is constant restricted to the leaves S x {i} of 
the time foliation. 

In itself the CMC gauge is only temporal, i.e. it fixes only the time foliation but 
does not fix the freedom by diffeomorphisms leaving the CMC foliation invariant. 

Assume that a vacuum cosmological solution (M, g) has a CMC Cauchy hy- 
persurface. Then, unless the solution is flat and of the form g = —dt^ -\- gp where 
gp is a flat metric over a three-manifold E, it must happen: a) there is a CMC 
foliation in at least a small neighborhood around it and b) if two CMC Cauchy 
hypersurfaces have different constant mean curvatures then they are disjoint and 
there is a CMC foliation inside the enclosed region on which k varies monotonically 
(for a proof see for instance [T7], [H], for a discussion see [25]). 

Let us explain how property a) is proved as that is relevant for the rest of the 
article. The basic tool is the technique of barriers |8J , it guarantees the existence of 
a CMC slice of constant mean curvature k between two slices with mean curvature 
functions fci and fc2 with sup ki < k < inf k2. This property settles the existence 
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of a CMC foliation around a CMC slice with fc 7^ 0. In fact starting at the given 
CMC slice S consider the foliation around it provided by the Gauss gauge with zero 
shift, i.e. with N = 1,X = 0. From equation ([7]) one gets k = > Thus k is 
strictly monotonic and we can conclude the existence of two barriers having mean 
curvatures ki and ^2 with sup ki < k < inf fc2- To show the existence of barriers 
around a state S with mean curvature zero proceed as follows. If J-^ \K\'^dvg ^ 
at the initial slice then we can solve the lapse equation — AA^+ \K\'^N = 1. In this 
case consider the normal field T — NT and consider the flow of geodesies starting 
at the slice and having velocities T there. This provides a gauge at least on a small 
time interval. The key fact is that at the initial slice we have A; = 1 and therefore 
k is strictly monotonic on a small time interval around the initial slice and in that 
gauge. This provides suitable barriers and therefore a CMC foliation around the 
initial CMC slice. If instead K = at the initial slice consider again the Gauss 
gauge starting from it. Let us restrict to the future direction. Unless K = over a 
small time interval (and in the Gauss gauge) we have from the equation k = \K\'^ 
that any other slice except the initial it is X 7^ and fc > 0. We can solve therefore 
the lapse equation at a time after the initial time and proceeding as in the previous 
case get a slice with inf k > 0. Doing the same in the opposite time direction we 
guarantee the existence of barriers and therefore the existence of a CMC foliation 
around the initial state. If it happens that if = on a small time interval (in the 
Gauss gauge) then from equation ^ we get that Ric = on a small time interval 
and therefore the solution is flat and of the form g — —dP + gp where gp is a flat 
metric over a three-manifold S. 

Fact a) above shows that any CMC state (with k ^ 0) induces a (at least 
short time) solution of the CMC flow, while b) shows that the CMC flow is unique 
and the maximal range of k forms a connected open interval over which k varies 
monotonically. The range of k depends on the topology of S and more in particular 
on the Yamabe invariant Y(Y,) of 10. It is conjectured [23] that if F (S) < then the 
range of k is (—00, 0) while if F(I]) > then it is (—00, -|-cxd). The CMC foliation 
is conjectured to cover the maximal globally hyperbolic solution if y(S) > while 
if Y{T,) < it is conjectured to cover the maximal globally hyperbolic solution 
towards the past but not necessarily towards the future. In this last case the CMC 
foliation should avoid the singularities but nothing else than that (every point lying 
outside the CMC foliation should have all its future inextendible time-like geodesies 
extinct in a uniform interval of proper time). 

Over a CMC foliation the function k is smooth except possibly at a maximal 
slice (i.e. a slice with fc = 0). The mean curvature fc is a smooth function at a 
maximal slice if \K\'^dvg 7^ at the slice. If we take the mean curvature fc as 

*The Yamabe invariant Y{T,) of a three-manifold S is defined as the supremum of the scalar 
curvatures of unit volume Yamabe metrics, where Yamabe metrics are metrics of constant scalar 
curvature minimizing the Yamabe functional. The Yamabe invariant is also called the sigma 
constant 1151 . [3|. 
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time (when k ^ 0) we get the lapse equation 



(8) - AN +\K\^N ^1. 

which is independent of the shift. 

Finally a note about the set of CMC solutions. Not every cosmological solution 
(Mjg) admits a Cauchy hypersuface of constant mean curvature [H] although 
there are sufficient conditions for that to be [S]. On the other hand every three- 
manifold admits so called Yamahe CMC initial states, i.e. states {gyjKy) with 
RgY = — 6, Ky = —5. Therefore there are CMC solutions in any S x M topology. 

A CMC solution over a manifold with non-positive Yamabe invariant and with 
the range of k containing an interval of the form (a, 0) is said a long time CMC 
solution (towards the future k exhausts its potential range). 

2.3 The CMC flow Equations in a general material setting. 

When matter is present, the Einstein equations 

Ric - ^Rg = SttGT, 

in flow form are 

(9) R~\K\^ + k^ ^leirGp, 

(10) VK = -SttGJ, 

(11) g = -2NK + Cxg, 

(12) K = -VViV + N{Ric +kK - 2KK) + CxK - 7V87rG(T - i<rgTg), 

where T(T, T) = p and T(r, .) = J and p is the average of the principal pressures 
(thus —p + ip = tr^T). The last term in the right hand side (RHS) of equation 
(fT2|) (involving T) should be restricted to S. The lapse equation for the CMC time 
t — k \s 

-AN + (47rG(p + 3p) -I- \K\^)N = 1. 

Some useful terminology. 

We have been using informally the notion of "control" . It is convenient to 
use this terminology indeed, and we will appeal to it many times in the rest of 
the article. In an informal manner, "A" "controls" "B" if B cannot degenerate 
when A is not degenerate. For instance, V and Qo control the "geometry" of 
CMC states. Quantitatively, (and being precise), a set of quantities Ai,...An 
control a quantity B if for every constant M there is N such that \B\ < N if 
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1^1 1 < . . . , \ An\ < M . Another more familiar way to mean control is through 
the expression \B\ < C{Ai, . . . , An). Still we will abuse slightly this definition an 
use it in a more flexible way. For instance we may say: 'V, V, Qo control r2", 
although properly speaking one must say: V and Qo control r2 + l/r2. The 

reader should find evident the meaning of each expression. 

2.4 Manifolds and Sobolev spaces. 

Given a subset Q of M" covered by a chart {x'^}, we denote by H^^y{Q) the 
s-Sobolev space of distributions with derivatives in until order s, and where the 
derivatives are the standard partial derivatives in the coordinate system {a;''}. 

We say that the pair is a -manifold if and only if A is an atlas 

covering E with transition functions in H^^^ (more precisely if {xa} is a chart 
from on a domain fj^ and X/j is a chart from a domain Vl' then a;J,(a;^) are functions 
on i?|^^j(ri n ^'))- We say that (E,^, g) is a W^^-Riemannian manifold if the 
entrances gij in any coordinate chart {xa} are in j(f^) where fl is the domain 
of the chart. 

Let {B(oa, r), a — 1, ... , m} be a covering of E such that every one of the balls 
B(pa,r) is inside one of the charts in A. Say {^Q,a = 1,...,to} is a partition 
of unity subordinate to the covering {i3(oQ,r)} (i.e. $q(p) = 1 if p G B{oa,r/2), 
^a{p) = if p e B{oa,rY, > and Ea^"(p) = 1 for all P e S). RecaU 
that in M" the Sobolev space is defined as those distributions for which the 
Fourier transform m(<j) satisfies that /u„(l + |'^|)^*|M(<?)p(i<; is finite. Thinking that 
the charts {xq,} are coordinates in (a subset of) R", consider the tensors U of rank 
(/, and having entrances in i?* for s <i and i > 2. Then the Sobolev space iJJ^ of 
tensors of rank (Z, V) is defined as the set of [/ = faf^a where Ua is a tensor 

as mentioned before. The Hj^ norm is defined by ||f/||//^ = X]q=i' llCa^alli?! ,• 
As topological spaces they are independent of the atlas A (as long as the atlas are 
compatible) . 

The H^-harmonic radius ri(p) at o in a i?*+^-Riemannian three-manifold (E, A, g), 
i > 2, is defined as the supremum of the radius r for which there is a coordinate 
chart {a;} covering B{o, r) and satisfying 




where in the sum above / is the multindex / = (ai, a2, as), and as usual /dx^ — 
{dx^)"^ (dx^)"'^ (dx^)"'''' . Both expressions above are invariant under the simultane- 
ous scaling g — X^g, 5^ — Xx^ and f = Xr. Observe that ii j > i > 2 then 

will always be an open set with smooth boundary. 
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'"j(o) < 'i'i(o). Define the H^-harmonic radius Ti to be the infimum of ri{o) for all o 
in E. The harmonic radius marks the scale at which the metric is seen in balls of 
radius one controlled in H]^^^^ by one. Atlas made of harmonic charts will be called 
harmonic atlas. Still one wants to have a more standardized notion of harmonic 
atlas. For that we define a H^- canonic harmonic atlas as one having the least num- 
ber of harmonic charts for which: every chart {xa} is defined on a ball B{o,ri/2) 
and can be extended to a chart on B{o,ri) having the properties (|13p - p^ . Also 
every ball B{o' , (3/8)ri) (for arbitrary o') must lie inside one of the charts. 

When we consider norms with respect to a canonic harmonic atlas A, it is 
necessary, in order for them to be useful, that the partition function with respect 
to which the norms are defined (see above) is not unnecessarily degenerate. More 
precisely the H'^^ynorm of a function with support in a chart {xa} in a H^- 
canonic harmonic atlas A does not have to be big if it can be small. One way to 
implement the criteria is to consider functions defined canonically as follows. 
Pick a C°° non-negative function x of one variable, being one inside the interval 
[0,3/4] and zero inside the interval [l,oo]. Given a chart {xa} in a _ff*-canonic 
harmonic atlas A, define Xa{x) — x(^|2^|)i where |a;| is the radial coordinate from 
the center Oq,. Finally define 



Let (S,(?) be a i/*+'^-Riemannian manifold, i > 0. As it is usual one can define 
the -ff-'-Sobolev norm of a tensor field U intrinsically using the covariant derivative 
V associated to g by ||f/||^, = /sEm=o iV^t/pdug. The H{ and the iJ^ Sobolev 
spaces are equivalent for i > 00 

All through the article we will consider pointwise products of tensors lying 
in different Sobolev spaces and whose product we would like to think in a third 
Sobolev space. All the mathematics we will need in this respect is the following 
fact (see [S] page 5 and references therein). Let U and V be two tensors over a 
iJ'+^-Riemannian three-manifold (S, g). Assume U is in and V is in , with 
i + 2>si>0, «-|-2>S2> cflll and one of the si or S2 is strictly greater than zero. 
Say s < min{si, S2, si + S2 ~ 3/2). Then 

(15) \\U*V\\H:,<C{r^.+2,Vol)\\U\\H^^\\V\\H^., 

where A is a, H'+^-canonic harmonic atlas. Thus for example if si > 2 we have 
independently of si, that the pointwise product J7 * is in H''_^ and its norm can 
be estimated from above by p^ . On the other hand if si = S2 = 1 the pointwise 
product of U and V lies in and its norm can be estimated by (jlSp . The reader 
should be aware these product properties are used intensively (specially inside the 
Section I3.3.2p , and most of the time without explicit mention. 



^''The case i = needs a bootstrapping argument to show that if V(/ = dU + T *U and U are 
in then dU is also in . 

^^The condition that < i + 2 for A; = 0, 1 is to guarantee that the Sobolev norms /f^* are 
well defined globally over E. 
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We will recur to elliptic estimates in a number of occasions. Except for standard 
elliptic estimates [16 and the elliptic estimates stated in Proposition [28] later in 
the article, all the rest will fall under the following proposition. 

Proposition 1 Say r < ri-)_2(o), i > 0, and let {x} be a harmonic coordinate 
chart covering B{o,ri^2{o)) and satisfying I113\) and We will consider elliptic 

estimates for first and second order differential equations separately. 

I. (Second order elliptic operators). Consider the differential equation 

(16) g'^d^djU'"" + + = /™. 

with A m W^^\{B{o,r,+2{o))), B m Hl,^^{B{o,r,+2{o))) and f m Hl^^{B{o,r,+2{o))) 
andwithj<i. Suppose in addition that \\A\\f^i+ii^gf^^ ^^^f^^^^ + \\B\\fji^ ^B{o,ri+2(p)) — 
C(r). Under these hypothesis we have that if U is a H^'^^ -strong solution of Hid]) 
then 

If B ^ we can get an estimate of one degree higher and JiTp is valid for j ~ i + 1 
as well. When i>2 the norm \\U\\hi on the RHS of equation ( [j7| ) can be replaced 
by the L^-norm of U . 

II. (First order elliptic operators). Consider the differential equation 

(18) G™a,C/" + A™C/" = /". 

Assume G^^diU^ is a first order elliptic operator with coefficients G™* in 

(i3(o, ri+2(o))) and ellipticity constants controlled by r. Assume also A is 
in 7J|^|(_B(o, ri+i(o))) and f is in _ff|^-j(i?(o, ri_|_2(o))) with j < i + 1 and finally 
that l!G||^j+2(3(^^^ (^^^^ + (o))) < C{r). Under these conditions 

we have that if U is a H^^^ -strong solution, then 

(19) ||C/||^. + i(B(„,./2))<^(^'^M)(l|t^llHj.,(B(o,.)) + 11/11^^^ 

If i > I then the norm \\U\\hi on the RHS of equation U9\) can be replaced with 
the L^-norm of U . 

Remark 1 The elliptic estimates above can be made global in a concrete way. In- 
deed let [Yi^A^g) be a i?*+'^-Riemanian three-manifold, where ^ is a 7f*+^-canonic 
harmonic atlas. Then, if we take Sobolev norms with respect to the atlas A we 
have, for instance the elliptic estimates 

(20) <G(n+2,G,l/oO(||C/||Hi, + 11/11^^). 

for second order elliptic equations of the form (fT6|) . The constant C bounds the 
norms of the coefficients A and B in H^^'^ and H''^^ respectively. Proving that 
the constant C in the elliptic estimate above can be made dependent only on G, 
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ri4-2 and Vol is not difficult. In fact and Vol control the number of charts in 
a canonical atlas and the iJ^^^j^-norms of ^q. It follows from the definition of 
given above that the global elliptic estimates are a direct consequence of the local 
elliptic estimates inside each one of the harmonic charts. 

Proof: 

I. Say i > 0. If we prove the case j = we can proceed by induction to prove 
the proposition for any j < i. Indeed assuming it is proved for j = jo < i — I, 
differentiate with respect to dk for k = 1,2, 3, use the result for j = jo and with 
dkU as U. Similarly the estimate when B ^ can be obtained from the estimate 
for j = i by differentiating the equation (flB)) with respect to dk, k = 1,2,3 and 
applying the result when j = i to the function dkU as U. 

We proceed with the proof when j = 0. Let us divide the proof according to 
the cases i > 2, i — 1 and i — 0. We will forget including the domains where 
norms are defined. When it corresponds, for example when using standard elliptic 
estimates, one must restrict domains to proper subsets. We will also forget about 
the subindex {x} in the norms. We write for instance instead of L'^^y Finally 
C(r) will denote a generic function depending on r. 

Case 1, i > 2: inspecting the coefficients A and B the result follows from the 
standard elliptic estimates ([IS], Theorem 9.11). In this case we can replace the 
77^-norm of U on the RHS of equation (|16p by its L^-norm. 

Case 2, i = 1: note that \\A*dU\\L2 < C{r)\\A\\H-2\\U\\m and \\B * U\\l2 < 
C{r)\\B\\ffi \ \U\\hi. Now apply standard elliptic estimates to the equation 

with the RHS as a non-homogeneous term and use the estimates before. This gives 
||C/||h^ <C(r,(7(r))(||C/||Hi 

as desired. 

Case 3, i = 0: this case follows by bootstrapping. First note the following three 
facts 

1. If U is in L" (a > 2) then Holder inequalities give 

2a 



\\B*U\\l, < \\B\\l2\\U\\l^, = 

\\A*dU\\Lff < \\A\\Le\\^U\\L^^, (3 



[2 + a)- 
6ai 



6 + ai 

2. Standard elliptic estimates applied to the equation 

with the RHS as a non-homogeneous term, give 

\\U\\h2^, < C{r)i\\B * Uh, + \\A * dU\\L, + + 
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3. Sobolev embeddings give [Hj: 

i) if 2/3' < 3 then < C{r)\\U\\H2.,' with a = 

ii) if 3 + /3 > 2/3 > 3 then \\U\\^,,,^_3 < C(r)||C/||^2.^, , 

iii) \\^U\\L^^ < C(r)||i7||^2., with ai = 

Start applying 1 and 2 with a = 6: from 1 we get (3 — 3/2 and with that /3 we get 
ai = 2, using those values in 2 we get 

\\U\\h2../2 < C{r){\\B\\L2\\U\\Le + \\A\\Le\\dU\\L2 + \\U\\ls/2 + II/IL3/2) 
< CirmB\\L2 + \\A\\Le)\\U\\H^ + \\U\\h^ + 

Now apply 3 i) with /3' = /3 - 1/6 = 4/3 to get a = 12 and 3 iii) with /3 = 3/2 
to get cki ~ 3. Returning to 1 and plugging in those values we get /3 = 12/7 for 
a = 12 and /3 = 2 for ai = 3. We use next 2 with the minimmn of those (3, i.e. 
/3 = 12/7 to get 

\\U\\H2..2fr < C{r)m\L2\\Uh^2 + \\AUe\\dUhr2,. + \\U\\h^ + \\fh2), 

and in turn 

\\U\\H2..2n < CirMBU^ + || A|Ue)||C/||^^..V3 + \\U\\h^ + 
giving thus 

||C/||h.,i2/. <C7(r,C•(r))(||C/||Hl + ||/|U0■ 
From 3 ii) we get \\U\\co.i/i2 < C(r)||J7||^2,i2/7. With that we get \\B * U\\l2 < 
C(r)||B||i2||C/||^2.i2/7. We use this in ^ with /3 = 2 to get 

\\U\\h^^^ < C(r)(||B|U2||C/||^2,i2/r + WAUeWdUUs + \\U\\h^ + \\fh2, 

giving 

||C/||ff2,2 <C(r, C(r))(|lC/|lffi +11/11^2), 

as desired. The elliptic estimate on first order elliptic operators follows exactly 
along the same lines as we did for second order elliptic operators. □ 

To illustrate the use of the estimates in Proposition [1] let us prove a well known 
fact that will be needed later. 

Proposition 2 At any point o on a compact H^^^ -Riemanninan three-manifold 
E (i > 0), n+2{o) and ||-R«c||^.+i^^(^_^^^^(^-|-|-| control 7-^+3(0) from below. As a 
consequence, for any compact H^~^^ -Riemannian three-manifold S, ||i?ic||^i+i^2) 
and ri^2 control r^+s from below. 

Proof: 

Pick a harmonic coordinate chart {x} covering i?(o, rj+2(o)) and having the 
properties (fT51) - P^ . In harmonic coordinates the Ricci tensor Ric has the expres- 
sion ^ 

Ricap --^g''^did-jgaf3 + d*g*d*g*g*g*. 
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The Elliptic estimates in Proposition [T] (where A — df,g*g*g* and _B = 0) show 
that ||5a/3||ff.+3(B(o,^,^,(o)/2))) is controlled by \\Ric\\Hi+\B(o^r,+2(o))) and ri+2{o), 
therefore ri+'i{o) is controlled from below by them too. □ 

Finally let us state the fundamental theorem of convergence in Riemannian 
geometry. Recall that the volume radius [2] at a point o is defined as 

^{o) ~ sup{r/vol{B{p,s))/s^ > fi, for all B{p,s) C B{o,r)}, 

where /x is an arbitrary small numeric constant. The volume radius of E is defined 
as the infimum of i>{o) for all o in E. 

Theorem 4 (Fundamental theorem of convergence (see and reff. therein)). 
Say (E, A, g) is a compact H^'^^ -riemannian three-manifold. Then the H'"'^^ -harmonic 
radius is controlled from below by Volg{T,), \\Ric\\Hi and Vg. Any sequence of met- 
rics with Volg{T,) and \\Ric\\iji uniformly bounded from above and Vg uniformly 
bounded from below has a subsequence converging (up to diffeomorphisms) to a 
-metric in the weak -topology. 

Remark 2 The use of the volume radius is not strictly necessary and it may well 
be substituted with other notions of "local flatness" (scaling as a distance) better 
adapted to use in General Relativity and more specifically in its dynamics. 

2.5 The Bel-Robinson energy and the space-time curvature. 

In this section we introduce Weyl fields and Bel-Robinson energies. We do so 
without further explanations or proofs. The reader can rely in [ISj for a detailed 
account. We will follow it in the presentation below. 

A Weyl field is a traceless (4, 0) space-time tensor field having the symmetries 
of the curvature tensor Rm. We will denote them by ^Vabcd or simply W. As an 
example, the Riemann tensor in a vacuum solution of the Einstein equations is a 
Weyl field that we will be denoting by Rm = Wq (we will use indistinctly either 
Rm or Wq). The covariant derivative of a Weyl field V^W for an arbitrary vector 
field X is also a Weyl field. In particular V^Wq = are Weyl fields. We will be 
using the Weyl fields W; with X = T, where T is the future pointing unit normal 
field to the CMC foUation. 

Given a Weyl tensor W define the left and right duals by *'Wabcd = ^fahinx 
and W*j,^^ = W^J^iez^cd respectively. It is *W = W* and *(*W) = -W. Define 
the current J and its dual J* by 

VWabcd = Jahc, 

When W is the Riemann tensor in a vacuum solution of the Einstein equations 
the currents J and J* are zero due to the Bianchi identities. This fact will be of 
fundamental importance latter. 
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The L^-norm with respect to the foliation wiU be defined through the Bel- 
Robinson tensor. Given a Weyl field W define the Bel-Robinson tensor by 

Qa6cd(w) = w,,„„w^™ + w:,_w%'7. 

The Bel- Robinson tensor is symmetric and traceless in all pair of indices and for any 
pair of timelike vectors Ti and the quantity Q{TiTiT2T2) is positive (provided 
W^O). 

The electric and magnetic components of W are defined as 

(21) Eab = WacbdT'T'', 

(22) Bab WacbdT'T''. 

E and B are symmetric, traceless and null in the T direction. It is also the case 
that W can be reconstructed from them (see [13], page 143). If W is the Riemann 
tensor in a vacuum solution we have 

(23) Eab = mCab + kKab " K a' K\, 

(24) ejBi, = Va/U, - "^bKac- 

The components of a Weyl field with respect to the CMC foliation are given by 
(i, J, k, I are spatial indices) 

(25) WyfeT = —^ij^'Bmk, *WyfeT = e^j"^Emk, 

(26) Wyfci — ^ijm^klnE"'" , *W,yfe; = f-ijm^klnB™"^ . 

We also have 

Q{TTTT) = \E\^ + 

Q.TTT - 2{EhB)„ 

Q,,TT = -{E X E),j -{Bx B\j + i(|£;|2 + \B\^)g,j. 

The operations x and A are provided explicitly later. The divergence of the Bel- 
Robinson tensor is 

V''Q{W)abcd =Wfc"S "J(W)„e„ + Wfe"2 "J(W)„rf„ 

+ * Wr/J*(W)„,„ +* Wr,"J*(W^)™e„. 

We have therefore 

V"Q(W)„TTT = 2E'\W)3{W\tj + 2B''^J*(W),Tj. 
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From that we get the Gauss equation which will be used several times during the 
article 

Q(W) = - / 2NE^\W)3{W)iTj + 2NB^^{W)r{W)iTj 

(27) Jt. 

+ 3NQabTT'n"''dVg. 

Hab = '^aTb is the deformation tensor and plays a fundamental role in the space- 
time tensor algebra. In components it is 

Ily = —Kij, TJiT = 0, 

LlTi — ^ , llTT — U. 

The next equations are essential when it comes to get elliptic estimates of Weyl 
fields. 

(28) divE{WX = {KA B(W))„ + JTaT(W), 

(29) divB{W), = -{K A E{W)), + r^.^iW), 

(30) curlBabiW) = E{VTW)ab + 1{E{W) x KU - ^kEabiW) + J„t6(W), 

(31) curlEabiW) = B{VTW)ab + ^(B(W) x K)ab - ^kBabiW) + J:t6(W). 
The operations A, x and the operators Div and Curl are defined through 

{A X BU = e^^\^U^eBdf + B)gab - ^{trA){trB)gab, 

{A A B)a = e^'A^'^Bac, 

{div A)a = VbA^^, 

{curl A),, = i(e„""V,^„6 + efe'^V^^^J. 

Finally wc mention a couple of formulas that will be used. If F a symmetric 
traceless (2, 0) tensor with 

{div V)a = VbVi = p, 

{curl V)ab = ^(ea'"VjKzb + e,'"^ViVrna) = <y, 

then 

(32) |Vy|2 + 3 < Ric, VoV> -\r\V\'' = £ \a\^ + \\p\\ 
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We have also 



(33) 



d^*d^iA) + 2div*div{A) = 2V*VA + TliA). 



where TZ has the expression Tl{A) = Ric o A + Ao Ric — 2Rm o A. The expression 
Rm o A is defined as {Rm o A)ab = Rmacbd.A'^'^ . 

The connection coefficients. 

Say T is the CMC foUation. Say X and Y are vector fields tangent to T and 
commuting with NT. In this setting we have 

K{X,Y) = - <VxT,Y >, 

VxY^ = Vx^*" - K{X,Y)T'', 



VxT^ = -X°-KJ', 



N ' 



where Pab = gab + TaT^ is the horizontal projector (it projects into J^), 

T 



and 



N 



2.6 Scaling and Cosmological scaling. 

Given a solution g of the Einstein equations, we say that A^g is the solution g 
at the scale of and we call A the scale factor. We say that a CMC state [g, K) is 
cosmologically scaled or normalized if fc = —3. This condition has a cosmological 
interpretation. We will see later that the Hubble parameter Ti can bo identified 
with =^ and in this sense a cosmological normalized state is one for which its 
Hubble parameter is equal to one. In general, space-time tensors, scale as for 
some weight s. The table below shows the scaling rules for some common tensors. 



g 


Ric 


T 


9 


K 


k 


N 


P 






A^g 


Ric 


T 


A^5 


\K 


k 
A 




A-V 




A-(2*+i)Q, 



3 The CMC flow, the volume and the BR-energy. 

Assumption: from now and for the rest of the article we will assume CMC states 
{g, K) have k and that CMC flows are defined on a range of k of the form (a, b) 
with —00 < a < 6 < 0. Also the future will be the direction in which the volume 
expands. 
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3.1 The reduced volume. 



In this section we assume that Y^E) < 0. The reduced volume that we wiU use 
was introduced and used systematically in [15] as a reduced Hamiltoniarl^ after 
passing to conformal variables[il. It is defined (up to a constant from [15] ) as 

Recall that in any flat cone the volume of the CMC slices grow as l/(— fc)'^. With 
that in mind we may interpret the reduced volume as a comparison between the 
volume of the particular solution and the volume of a flat cone at a given A0. We 
enumerate below a list of properties of V [15] . 
1. The derivative of V with respect to k is 



(34) ^ = -ih' I N\k\-^dv, = -ilr I (1 - ^)dv. 



dk 3 J Y; 3 J Yi 



^9- 



When V = then N = ^,V" = and V" = -fc^/g N^\Ric\^dvg. 

2. From the property before we get that V is strictly monotonically decreasing 
along the future CMC flow unless is constant in which case the solution is a 
flat cone. 

3. V is scale invariant. 

4. The inflmum Vmj of V in the phase space of all CMC states is given by 

V„/ = inJ{V{g,K)/{g,K) is CMC} = (-iy(S))i. 



where 1^(S) is the Yamabe invariant of S. 
Under zero shift we have the remarkable fact that 

= I NR + Nk^dvg = —VoL = -^V 
< 'MV, 



^^In |15l the reduced volume is named the reduced Hamiltonian. 

^^In the context of the long time a quantity similar to the reduced volume was also used in 
[l]. More precisely, the "reduced volume" in [l] was defined dynamically: say {{g, K), (N, X)) is 
a solution to the CMC flow and tj. = distCSf^^ , Sj.) where dist is the Lorentzian distance between 
the initial CMC Cauchy hypersurface Sj; and Sj. , then define V = ^ °'g^^^ . It is proved [T] that 

it is a strictly monotonic quantity unless is constant in which case the solution is a flat cone. The 
reduced volume here defined has the advantage that it is defined on the set of all CMC states 
whereas the one in [A] is defined along solutions. 

^■'The flat cones are the solutions with maximal rate of volume expansion. 
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where TV = -^^^^ is the lapse associated to the cosmological time t = — 3/fc and 
g — dg/dt. To obtain the first inequahty we have used the fact that iV < 1 and to 
get the equality after the fourth term we have used the lapse equation ([5]) multiplied 
by As V is monotonic along the future direction in a CMC solution of the 

Einstein vacuum equations, the RHS of equation is bounded above by 3tV(io) 
where to is some initial time. As for the traceless part of the time derivative of g 
we have 

2 



implying 

1 Z"*! \\9\\l2 



4 L. t^ 



dt < V(<o) - V{ti). 



3.2 The Priedman-Lemaitre equations and a cosmological in- 
terpretation of the reduced volume monotonicity. 

The Robertson- Walker cosmologies implement the cosmological principle in its per- 
fect form. As such it characterizes the universe by a radius a, matter density p and 
the pressure p. The space-time is modeled by a metric 

-dt^ + a^{t)gK:, 

on a manifold E x M = M'^ x M and where gjc is a metric of constant sectional curva- 
ture equal to /C. It follows that the stress energy-momentum tensor is of the form 
Tah = {p+p)TaTh+pgab with p and p depending only on time. The compact FLRW 
models are simply space-like compactification of the models above. In the IC — 
case compact Cauchy hypersurfaces obtained are flat while in the /C = — 1 case 
are hyperbolic. The Einstein equations are equivalent to the Friedman-Lemaitre 
equations that we present below in integral form to make a closer connection with 
the Friedman-Lemaitre equations for arbitrary solutions that we describe later 
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2 _ /g(167rGp)dt;g _r(}^\^ 
a" - J^{4ttG{p + 3p)dvg 



a 3V 

In the formulas above Vq is the volume of the manifold S when it is given the 
unique hyperbolic metric or some flat metric according to the type of the model. 
In a general cosmological solution (not perfectly homogeneous and isotropic) the 
cosmological parameters are defined in volume average [22] , [TOlf^ . For instance the 
radius a is defined as 



""Although the averages in |22| and |10| both average in volume, they are not precisely the 
same. 
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and the proper time by 

dr /j. Ndvg 



dk V 
The Hubble parameter Ti, is computed as 

■N - - - —-V--( f Nhrl N ^ _ -fc 

This inspires to give the name of "Hubble-gauge" to the CMC gauge, as observers 
at the same "instants of time" would measure the same Hubble parameter H (i.e. 
it is constant along leaves of the CMC foliation). Define p = T(T, T) and — p+3p = 
Tabg"'' (note that p is the average of the principal pressures). Then the Friedman- 
Lemaitre equations in this general setting are 

J^Rdvg ^ J^{16nGp+\k\^)dvg 



6V 6V 
J^MiiirGip + 3p) + \k\^)dvg 



(36) H=-- 

Where Af — ^ (bar denotes volume-average) and has average equal to one. All 
the derivatives above are with respect to the averaged proper time r. The first FL 
equation (equation (|36|) ) is just the average of the energy constraint 

16np = R- \k\^ + -k^. 

3 

To get the second FL equation (equation (|57|) ) we observe that 

(38) (^)'^^_(^)2 = ^_^2^ 

a a a a 

and 

(39) C-y ^dn^^idj^^ V 



a dr 3 dr 3 /j^ Ndvg 
On the other hand from the lapse equation we get after integration 

(40) / N{4:nG{p + 3p) + \k\^)dvg -3n^ I Ndvg. 

Equations (|38p. ((39| and (|40|) together give the equation p7)) . Up to a constant 
the reduced volume is V = (3ai7)'^ and we see that 



{V^y ^ {3aH)' = 3a" = - 



a^^N{A^G{p + 3p) + \k\^)dvg 



V 

If there is no matter present, i.e. p = p — the previous equation is directly equiv- 
alent to the monotonicity formula for the reduced volume. This implies that 
the monotonicity of the reduced volume and the well known universe deceleration 
are equivalent properties. Note that the reduced volume is monotonic in presence 
of matter at least if p -|- 3p > 0. That is the case if the average of the principal 
pressures is positive or if the mass density dominates over the (possible negative) 
pressure. 
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3.3 Controlling the states {g, K) at a given time. 

This section is devoted to study the intrinsic quantities V, Ej = X]i=o Qi^i) 
and p (if matter satisfying the dominant energy condition is present) as variables 
controUing the states {g, K) at a given time t = k. 

Assumptions. All the results presented in this Section fSection 13. 3p . except 
for some exceptions explicitly indicated, will be stated in the context of vacuum 
solutions of the Einstein equations. Also, all the analysis will be for general three- 
manifolds S, except when we indicate that a result is restricted to manifolds with 
F(E) < 0. We will assume that all quantities (for instance g, K or W^) come from 
a C°° solution g. In this way we avoid justifying certain operations on tensors 
that would require justification in the case of low regularity. Still the results of 
this section extend to the natural regularity of each statement or proposition. For 
instance Proposition [3] is valid for states {g, K) with x regularity. The same 
holds in Proposition [TH where to exemplify, a footnote was added on how to prove 
the result for states with x regularity. All the time C will represent a generic 
constant (the same C can represent different constants in different lines). If the 
constant C is numeric it will be stated explicitly. 

The main result in vacuum will be the following. 

Lemma 2 (Sobolev norms vs. Bel-Robinson norms) Say I > and say T, is a 
compact three-manifold. Then the functional (defined over cosmological normalized 
states {g, K) ) 

\\{g,K)\\Bn = ^+V + £j, 
controls the H^'^^ -harmonic radius rj^2 ^'^'^ ^1^^ H^'^^ -norm of K . 

The proof of this lemma (at the end of this section) is conceptually divided into 
several propositions many of them having however independent value. As it turns 
out there are a number of intrinsic estimates, namely estimates involving only 
intrinsic norms i?*, worth to be mentioned separately. Intrinsic estimates give more 
information than what elliptic estimates give. The elliptic estimates are stated in 
terms of the norms Hj^ where A are harmonic atlas. To highlight the difference 
between them, we break the section into two subsections. 

3.3.1 Intrinsic estimates. 

We start studying the control of the first order Bel-Robinson energy Qq over 
states {g,K). 

Proposition 3 Say Yi is a compact three-manifold. Then Qo (md jfcpHA'll^g con- 
trol \\yK\W2, 11-^11^4 and \\Ric\W2. More in particular we have 

(41) (/ 2\Vk\'' + \k\Uv,)-- <C{\k\\\k\\L.+Ql), 

Jm 
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where C is a numeric constant. 

Observe the absence of the volume in equation (|4T|) and that aU norms involved 
are intrinsic. 

Proof: 

Substituting Ric = E-kK + KoK,K^K + ^g and V ^ K in equation ^ 
we get 

/ \Vk\^ + ^\K\'^-k < k,KoK > -^\k\'^ + 3 < E,kok > dvg = [ \Bfdvg. 
This equation gives the bound 

(42) / \Vk\' + \k\Uvg <C f m^\K\^ + \k\\K\' + \K\^\E\ + mdvg, 
je Je 

Observe now that the inequalities 

/ \k\^{\E\' + \B\^)Uvg<i[ \k\Uvg)^Ql, 

JS JE 

/ \k\'dvg<{f \k\^dvg)H [ \k\^dvg)K 

i/ 5j •-'2 £ 

transform equation (j42p into 

2\\vk\\l. + lli^lli. - C(|fc|||i^|U. + Q|)ll^lli^ - C{\k\'\\k\\l, + Qo) < 0. 

Now make ^ 2 \VK\'^ + \K\'^dvg, a = (|fc|||ii:||i2 + Q'^) in the last equation. 
We get 

- Cax - Ca^ < 0. 

Solving for x in the inequality above we get equation (|4ip which finishes the proof. 
□ 

In presence of matter satisfying the dominant energy condition (see [18) for a 
discussion) we have the following modification of the last proposition. 

Proposition 4 Say Y. is a compact three-manifold and say (g, K) is a cosmological 
normalized state over a Cauchy-slice on a solution (M, g) of the Einstein equations 
and with matter satisfying the dominant energy condition. Then ||_ftr||^2, Qq o,nd 
||Gp||l2 control \\\JK\\]^2, \\K\\]^i and \\Ric\\i^2. More precisely we have 

(43) {jj\Vk\^ + \k\UvgY^ <C{\\k\\L-^+Ql + \\Gp\\Li) 
where C is a numeric constant. 
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Proof: 

If the energy-momentum tensor of matter satisfies the dominant energy condi- 
tion then Ttt > |Tij| for any i,j — 1,2,3. Recall that the space-time curvature 
Rm is decomposed in terms of Wq and Ric (with an expression linear on each) and 
that Ric — &ttG{T — ^tr ^Tg). We have then that the norm of the components 
of Rm are controlled by Qo + HGplUa- Using these facts the proof goes exactly 
parallel to the one of Proposition [3l but this time using instead the identities 

RmiTjT = Ricij - KimK"^ - kKij, 
RmmTij = (d^ K)m^j■ 
R^16nGp+\k\^ - ^fc^ 
VK = -SttGJ. 

□ 

The next proposition relates ||i^||^2 with V — Vm/ or V depending on the sig- 
nature of the Yamabe invariant 

Proposition 5 . Say Yi is a compact three- manifold. Then 

i) if > 0, \k\Hv, < c\k\'^vHj^ \k\^dv,)i, 

ii) if Y(S]) = 0, \k\^ \k\^dv, < C\k\Hv - V.nf)Hk \K\'dv,)i, 

iii) if Y(I]) < 0, \k\^ \k\^dvg < C{\k\{V - V,nf)+ 

+\k\Hv~v,^f)Hj^\k\^dv,)^). 

where C is a numeric constant. 



Proof: 

i) and ii) {Y{Y^) > or Y{Y,) = 0). This is immediate from the formula 

\k\' [ \k\^dv,<\k\H\k\'voiiY))H[ \k\Uv,f^. 

iii) y(S) < 0. Assume k — —3 and let gy be the unique Yamabe metric of 
constant scalar curvature Ry = —6 in the conformal class of g. If <? = (f'^^gy then 
(j) is determined by 

(44) - A,,^+^^~^r'\Kfy + ^k'cb' - 0, 

where A = V^. The maximum principle implies (putting the values Ry — —6 and 
k = —3) that 

6(0L„-0mm)>0„i|if|^>O, 

which makes (j> > I. Then observe that 

-y(S)<-i?y(/ ldvy)i, 
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where dvy = dyg^ ■ This gives 

0<6^(/ (/)^ - 1 dwy) < 6^ / (jj^dvY -{-Y{Y.)Y- 

= (^fc2Fo?(i])i)i-(-y(i]))i. 

Therefore 

\fdvY < C(V - V„/), 
for A: = 1, . . . , 6. Integrating equation pi)) . we get 



(45) G {(/)- (f)dvY = I <?!>"1-ftr|f rfwy. 

Observe that 



(t)'^\K\i^ dvy = / (t)(l3-'^\K\i. dvY 

/ (j^-^'ml dvy + / (0 - dvy, 

je je 



and 

(46) / (0 - l)r'l^l|- rft'y = / (0 - dvy 
je Je 

< ( / (</.- 1)V / r'"l^|y dvy)'^. 

Je Je 

On the other hand note that 

(47) I / " 1)'(/)^ dwrl < / - 1| + 2|05 _ i| + _ ^1^^^ < ^(-^ _ 

JE JE 

Putting together equations (|45l) . (|46|) and (|47|) we get 

(48) <c((v-v„/) + (v-v™/)3||x|i2,), 

which after scahng finishes the proof. 

Combining Propositions [3] and [5] we get 

Proposition 6 Say is a compact three-manifold. Then ifY{Y,) > Q we ha' 

2\VK\^ + \K\Uvg < C{\k\V + Qa), 
while z/y(E) < we have 

^ 2|Vif |2 + \K\^dvg < C{\k\{V - V,nf) + Qo), 
where C is a numeric constant. 
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We also get 

Proposition 7 Say is a compact three-manifold. IfY{T,) > we have 

I \k\''\k\''dvg<C{\k\V + {\k\VQ^Y^), 

while ifY{E) < (same for = than for < Oj 

/ \k\^\K\^dvg < C{\k\{V - V™/) + (|fc|(V - V„./)Qo)^). 

where C is a numeric constant. 
Proof: 

Combine equations in Proposition [5] and equation ((4T|) . Making x — |fc|||if||i2 
and a = |fc|^(V - H„/)^ if < or a = \k\^V if Y{T.) > wc arrive at the 

inequality x"^ — Cax — Ca^ — CoQq < 0. From it we get x'^ < C{a'^ + oQq). □ 

A direct consequence of the propositions above is 

Proposition 8 Say is a compact three-manifold, then V, \k\ and Qq control 
\\Ric\\i^2. In particular we have 

\\Ric\\l.<C{\k\V + Qo), 
where C is a numeric constant. 
Proof: 

Use Ric ^ E - ^^k + K o K - ^\K\'^g together with the Propositions Eland Ha 
Using the energy constraint R — \K\^ — ^k^ and Proposition [H] we get 

Proposition 9 Let T, be a compact three- manifold. Then, V, \k\ and Qq control 
the scalar curvature in the following way 

( \WR\^ +R^dvg<C{\k\V + Qo), 

where C is a numeric constant. 

Note that |Vi?|3 and R^ scale as a distance^'^ . 
Proof: 

Squaring the energy constraint and integrating we obtain 

j R^dvg < J \k\^ + ^k^dvg < C{\k\V + Qo). 

where in the last inequality we have used Proposition [6] On the other hand, 
differentiating the energy constraint we have |Vi?|^ < C|V-ftr|3|-ftr|^. Integrating 
and applying the Holder inequality we obtain 

/ \yR\Uvg<c{[ \vk\^dvg)l{[ \k\Uvg)K 

Jy. Jt, 
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and if we apply Proposition [H] over each one of the factors on the RHS of the last 
equation we obtain 

/ \VR\Uvg < C(|fc|V + Qo), 

as desired. □ 

Because of the monotonicity of the reduced volume in the future direction we 
have 

Proposition 10 SayY{Y.) < 0. Then along the future CMC evolution Qq controls 
WmcUp \\K\\hi and ||i^|U4. 

In the presence of matter satisfying the dominant energy condition we have 

Proposition 11 Say Y. is a compact three-manifold. Then (over cosmological nor- 
malized states) Qq, \\Gp\\]^2 andV — Vinf control \\K\\fji , \\K\\]^4 and \\Ric\\i^2 . In 
particular ifY{T,) < we have the formula (same for = than for < 0) 



Is 

while > we have 



\K\''dvg < CiiV - V,nf) + (V - V,nf)HQS + \\Gp\\l^)) 



\K\'dvg < Civ + VHQS + WGpWli)), 

where C is a numeric constant. 
Proof: 

y(E) < 0: the proof proceeds in parallel to the proof of Proposition [5] but this 
time instead using the equation 

and then the equation (|43p . 

> 0: use equation (l43l) in the formula 

/ \KfdVg<C\k\^VH[ \K\^dVg)K 



and after making x = ||A'||i2 solve for x. □ 

The lapse has important natural properties that we describe in the proposition 
below. 

Proposition 12 Say T, is a compact three-manifold. Then, ||iV||ioo < and 
\\N\\f}i is controlled by l/\k\ and V. In particular we have the bound 

(49) ^ I VTVp + \K\'N' + ^N^dvg = J^N< ^V. 

Proof: 

From the maximum principle we have N < 3/\k\'^. To get equation (|49p multiply 
the lapse equation by N, integrate, and use the estimate for ||Af||Loo. □ 
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3.3.2 Estimates of elliptic type. 



Assumptions in Section \3.3.2\ All through this section we will assume {g, K) is 
a cosmologically normalized state. 

We turn now our attention to the study of the influence of the higher order 
Bel-Robinson energies over states. 

We need a notion of Sobolev norm for space-time tensors and with respect to 
the CMC foliation. We do that in the following way. Consider a space-time tensor 
^ai,...,a, of rank (l, 0). For any subsequence / = (ii, . . . , i„) (n < I) of the sequence 
(1, . . . ,1), and the obvious compliment subsequence /, define 

Ti = Ta,^ . . . Ta,^ 

and 

where P^^^ = + TaT'^ is the horizontal projector. We can decompose U as 
U= J2 (^/<U,T,>)(T,), 

n— 1,...,/; \I\—n 

where < U,T/ > is the contraction of U and Tj. For each summand, the factors 
on the left are horizontal, while the factors on the right are vertical. For instance, 
for JJab we have 

Uab = P/Pb^'Va-b' - P/{\Ja'b'T''')n - Pb'' {U a' b'T""' )Ta + [U a'b'T'^' T'' )TaTi,. 

The same decomposition holds for tensors of arbitrary rank {1,1'). Now the H^- 
norm of U on a slice S of a CMC foliation is defined as the sum of the _ff^-norms 
of the tensors Pj < XJ,Tj >. We will be using this convention somehow implicitly 
all through and without further comments. 

From now on — V^Wq where T is the future pointing unit normal to the 
CMC fohation. 

During the proof of the propositions until the end of this section, we will use the 
notation iJ* instead of H'^^^y{B{o,r)) which is the one used inside the statements. 

Proposition 13 Say T, is a compact three-manifold. Then, (the data) v, V and 
Qo control ||iV||^2^ where A is a H^-canonic harmonic atlas. 

Proof: 

By Proposition [T21 is controlled by V and Qo- By Proposition [51 \K\'^ 

and therefore \K\'^ are controlled in by the data. The result then follows by 
Proposition [1] (I). □ 
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Remark 3 One can get an estimate for the intrinsic norm \\N\\h-2 in terms only 
of V and Qo ( i. e. without involving the volume radius ) if the Ricci curvature is 
hounded below. Indeed that follows from 

[ \VVN\^ + Ric{WN,VN)dvg = [ {ANfdvg, 

the use of the lapse equation in the RHS of the equation above and finally Proposi- 
tionsU^ and\^ 

The proposition below shows that for cosmologicaUy normahzed states, i', V 
and Qo control the norm of where A is a _ff^-canonic harmonic atlas. 

This implies in particular that the infimum of the lapse is never zero even for states 
with low regularitjlll. As a corollary, we get that ly, V and Qo control the H\ norm 
of the deformation tensor 11. 



Proposition 14 Let T, be a compact three-manifold. Then v, V and Qq control 
||1/A^||^2 where A is a H'^ -canonic harmonic atlas. In particular they control 

Proof: 

Multiplying the lapse equation by 1 /N'^ and integrating gives 
ViVP 1 . f \K\^ . , f . f 1 



This shows in particular that ||A^||/fi is controlled by V and Qo. We multiply now 
the lapse equation by and integrate, it gives 

By the Sobolev embeddinj^ Il\ ^ L\, the RHS is controlled by v, V and Qq. 
We will use this estimate below. Consider the Laplacian of 1/iV. We compute 



N N N ' ' 

We have then the elliptic non- homogeneous equation for 1 /N 

1 1 |V7V|M 



^'^We will operate assuming a priori regularity of 1/A'^. This is indeed guaranteed from the 
assumption made at the beginning of Section 13.31 To show that in this case the estimate also 
descend for states {g, K) with X regularity proceed as follows. Smooth out the coefficient 
l-R'P to get, from the maximum principle applied to the lapse equation, an upper bound on 1/N . 
This gives the necessary regularly to operate. As the estimate on the L°° norm of 1/7V that one 
obtains (having smoothed depends only on u, V and Qq it passes to the limit when the 

smoothing is undone. 

^^It is important here that the norm of the embedding is controlled from above by u, V and Qo- 
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where, to the effect of applying eUiptic estimates, we are thinking \\/N\'^ /N as a 
factor in front of the variable We know that \K\'^ is controlled in L^. The 

result then follows by Proposition [1] (I), if we show that \VN\'^/N is controlled in 
L^. We compute 

As was shown above, the first factor in the RHS of the previous equation is con- 
trolled by I/, V and Qq. The second factor is controlled by i^, V and Qq by the 
embedding ^ and Proposition [121 ^ 

Proposition 15 Say I > 0. Let f < r < r 1^2(0), and say {x} is a harmonic 
coordinate system covering _B(o, r/+2(o)) cind satisfying il3\) - (l^\ l. Then (the data) 
\\^o\\Hl^^(B{o,r)), \\^i\\Hl^^iB(o,r)), \\K\\ ^(^B{o,r)) , r and r control \\ Ric\\ j^i+i ^^g^^ -^^ 
and \\K\\^i+2f^gf^^ -^y In particular, they control r 1^3(0) from below. 

Proof: 

The proof proceeds studying the equation ([55)1 (making K = A) to get the 
estimate on K and an appropriate elliptic system f ((53| - (|54p ) to get the estimate 
on Ric. From equation (|24p we have 

and therefore 

(50) d^*d^i^ = -e/V,B,™ - eJ:V,Bij = -2curl{B). 
From equation pop we have 

(51) curl{B)=E{\/TW)+'^{E^ K)-]^kE. 
Equations ([50)) . dHU) and ((551) give the elliptic equation 

(52) 2V*VX = -nik) - 2(£;(Wi) + ^{E xK)^ ^kE). 

Recall Ti-iK) is a linear expression in K with coefficients involving only Ric. We 
consider the case / = first. This case in turn is the only one that demands a 
special treatment. In order to apply Proposition [1] (I) to the elliptic equation ([5^ 
we need first to obtain control on the H^-j^^y^ymycva. of K. This estimate follows 
from standard elliptic estimates on the elliptic system 

v^k,, = 0. 
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We apply then Proposition [T] (I), on the equation to get the control on the 
i??-norm of K. To get the estimate on the i?^-norm of Ric we get first an estimate 
for the H}-, ,„-norni of E. From it and the equation E = Ric + kK — K o K 
the estimate on Ric follows. To get the estimate on E apply standarqifl elliptic 
estimates on the elliptic system 

(53) curlE = -B{Wi)-'^{BxK) + ^kB, 

(54) divE^{KAB). 

Now we treat the cases / > 0. We note that as the harmonic chart {x} satisfies 
(fT3|) - (IT4|) we have control on |j_Ric|j//i by the data, and as was mentioned before 
the coefficients of K in the expression Tl{K) involve (linearly) only Ric. Thus we 
can apply Proposition [T] (I) to the elliptic equation ([5^ to get the desired control 
on ||K||^/+2. Similarly applying Proposition [1] (II) to the elliptic system ([5^- 
we get control on ||i^|| ^^z+i , and therefore on ||i?ic||j:^7+i from the equation 
E = Ric+kK -KoK. ^ " □ 

Proposition 16 Say I > 0. Let f < r < r/4.2(o), and say {x} is a harmonic 
coordinate system covering _B(o, r/+2(o)) and satisfying U3\} - (T4\ l. Then (the data) 
l|Wo||H|^j(s(o^r)), |lWi||^7^^(B(„ ,,)), \\k\\L2(^B{o,r)), \\N\\^(Bio,r)), r and r control 

W^W H'+fiBio.f))- 

Proof: 

We consider the case 1 = 0, the cases / > 1 easily follow by induction. By 
Proposition [T5l ||i?ic||//i^^ is controlled by the data. Therefore ||.g||H3^ is 

controlled by the data. By Proposition [15] too, is controlled in Hff^r)/2- 
can apply then standard elliptic estimates on the lapse equation 

-AiV+ \K\'^N = 1, 

to get that I lA^ 1 1/^3 is controlled by the data, as desired. □ 

Proposition 17 Say f < r < r/ 4-1(0) with I > I and say {x} is a harmonic 
coordinate system satisfying l7^) -|j^[j and covering i3(o, r/+i(o)). Then for any 
{i,j) satisfying < j < / and 1 < i < I , ||Wj||^i ^(B(o,r)) controlled by (the 

data) \\^j + l\\Hl-liB(o,r))' l|W,l|^-|(B(o,.))' l|Wo||ff-J(B(o,r)). 1 1 ^1 1 1 ^j- 1 , 

l|J(Wj)|lj^.-i(B(o^^)), r, r and \\k\\L2^B{o,r))- 
Proof: 

^*To apply standard elliptic estimates we note that as was proved before the coefficients involving 
K are controlled in , , , . . 

r+3{r — r) /4 
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Think the elhptic system (P5)l - (l3ip as a first order elhptic system of the form 
(fT8)) with U = {E{Wj),B(Wj)). By Proposition [15] the coefficients A™ which 
involve only K are controlled in i/|^|(i?(o, (f + r)/2)) by the data. The result then 
follows by applying Proposition [T] (11) to the elliptic system ([28|) - ([29l) . □ 

Remark 4 Applying Proposition [171 when j = we deduce that ||Wo||jyi is con- 
trolled by ||Wo||^i-i and ||Wi||^i-i, f, r and ||i^||i2. This tells essentially that 
one can replace the data ||Wo||j:^i-i, ||Wi|j by the data ||Wo||^i, ||Wi||^i-i 
inside those statements whose hypothesis contain the data r and ||_ft'||i2. This 
Remark will be used later (sometimes implicitly). 

We prove next, in Proposition 1 181 an inductive formula for the currents J(Wi) and 
then in Proposition [191 estimates on the time derivative of the deformation tensor 

n. 

Proposition 18 Say j > 1. J(Wj) has an expansion of the form 

(55) j(Wj) = ^(v™^n)"i * • • • * (v^'-n)"' * n' * vw^ 

(56) +^(V™in)"i *---*(V™=n)"^^ *n'* V|,(T*Rm*W^), 

where every asterisk * is some tensor product and each expression of the form 
(V™n)" is a ^-product o/ VtII with itself n-times. The sum on the RHS of Ii55\) 
is among sequences ((mi, tii), . . . , (mg, n^), I, k) with k < j ~1, nii > 1, . . . , > 1 
and "-^ (1 + fTT-j) + I + k ~ j , while the sum 156]) is among sequences 
((mi, rii), . . . , {ms,ns),l, q, k) with mi > 1, . . . , m-s > 1 and Yl,-j nj{l + fhj) + k + 
l + q = j-l. 

Proof: 

First note that 

(57 

= n * VWj + T * Rm * + VtJ(Wj). 

Now for j = 1 we have J(Wi) — 11* VWo+r*Rm*Wo which agrees with the form 
of the formula above ( ([55]) . ([5S]) ) . For j > 1 we proceed by induction. Assume the 
J(Wj) has the desired expansion. The first two terms in the RHS of equation ([57| 
are of the desired form, so it remains to prove that VTJ(Wj) is of the desired form. 
Terms in ([55]) and ([55]) are characterized by sequences ((mi, rti), . . . , (m^, Ug), I, k) 
with nj{l + rrij) +1 + k = j and ((mi, n-i), . . . , (m^, hs),l, g, k) with J^j + 
rhj ) + k + l + q = j — 1. We show that the Vt derivative of any term of the form ([55]) 
or ()56p gives terms of the same form characterized by sequences adding j + 1 for 
terms of the form ([55]) and j for terms of the form ([5S]) . Let us consider derivatives 
of terms of the form ([55]). The derivatives VT(V"n)" = n(VJ'+^n) * (V™n)"-i 
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transform the pairs (m, n) (inside a sequence characterizing a term of the form [55]) 
into pairs (m + 1,1), (to, n — 1) but leaving the rest of the sequence unaltered. As 
1(1 + 1 + to) + (jt. — 1)(1 + to) — 1 + n(l + to) the new sequence adds j + 1. Similarly 
the derivatives VtH' — ZVtII * II'^^ transform the I (inside a sequence) into the 
pair (1, 1) and the number I — 1 but leaving the rest of the sequence unaltered. 
Again in this case as 1(1 + 1)+^ — 1 = ^ + 1, the new sequence adds j + 1. Finally 
the derivatives VtVW^ are 

(58) VrVWfe^ VWfe+i+n*Wfe + r*Rm*Wfe, 

which transform the number k (inside a sequence) into three new sequences. One 
with the new number fc + 1 and leaving the rest of the sequence unaltered, thus 
adding j + A second with the munbers k and I = 1 and leaving the rest unaltered, 
thus adding j + Finally a third of the kind ((56)) which has values I = I, q = 0, 
k = k and {{mi, hi), . . . , {rhs, hs) = ((toi, ni), . . . , (tos, ris)) adding j as desired. 
The analysis of the Vt derivatives for terms of the form (j56p proceeds exactly in 
the same fashion. □ 

Proposition 19 Say f < r < r/+2(o); ^ > 0. Let {x} be a harmonic coordinate 
system satisfying the conditions il3\) - [T^ and covering _B(o, r/+2(o)). Then, for 
any {m,i) satisfying, to > 0, i > and 2<TO + i</ + 2, ||V™n||^i (B(o,r)) 
controlled by (the data) ||Wfc||^(m-fc)+i-i for k — 0, . . . ,m, 

{x} \ \ 1 J / 

W^T^W Hi+l"'-''>+\B(o,r))' for k = Q,...,m, \\k\\L2(^B{o,r)), r and r. 
Proof: 

The proof proceeds by induction. We first observe that the cases comprising 
those (m, i) such that m — Q and 2 < i < I + 2 are proved by Propositions 
1151 and 1161 The induction process will be as follows: assume the proposition is 
proved for all (to, i) with 2<m + i<I + 2 and m < fh. This assumption will 
be referred in what follows as the inductive hypothesis (IH). Under the inductive 
hypothesis we will prove that the proposition is valid for to = to + 1 and all i with 
2<i + m+l<I + 2. In this way we cover all (to, i) with 2 < m + i < I + 2. From 
now on we assume (to, i) =^ {fh + l,i) with 2 < i + fh + I < I + 2 and the data for 
(to, i) = {fh + l,i). 

Observation 1: A crucial observation which is easily checked is the following: 
the hypothesis of the Proposition for to = to + 1 and i = i, contain the hypothesis 
of the Proposition for all (to, i) with m < fh and 2<i + m<fh + i + l. In this 
way we have that || V^n||^i+(™_ii)+i , with < (5 < to, is controlled by the data of 
the Proposition for (to, i) — {fh + l,i). 

The following commutation relation will be used recursively. 

(59) VT{VaUb) = Va{VTUb) + T^U^Rm^aM ~ n^'V.Ub. 
From it, we get 

(60) V^+^Uab = V??(Va(VTTh) + T-T^'Rm^abd - U^'U.b). 
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We treat each one of the three terms that appear on the RHS of the last equation 
separately. We treat first the last term. We make the expansion 

From Observation 1, we get that at each summand, V^II and V^II, are controlled 
in and respectively. We get therefore that the full expres- 



sion is controlled in Hl^. We treat next the second kind of terms in equation ((60|) . 
We make the expansion 

We will treat this term using the following Fact, which, as the other Facts to be 
stated later, are going to be proved after the main argument is finished. 

Fact 1 : terms of the form for 1 < 6 < fh + 1 are controlled in iJp^'™ "^^^^ 

by the data. 

We know that ||W^||^m+i-^ is controlled by the data. It follows from this and Fact 

[T] that the second kind of terms in equation ([60|) are also controlled in H^. We 
discuss now the first kind of terms in equation (j60p . namely the terms of the form 

V??(Va(VTTb)). 

We would like to pass the Vt's on the left of this expression to the right of Va- 
We will show that every time a Vt is moved past of Va we generate a pair of 
terms that are seen to be controlled in H^. We write 

V??(V„(VtT,)) =V™-l(Va(VT(VTn)) + T^(VTT'^)RmeaM 
-n,^(Ve(VTT,))). 

We state now Fact [H and Fact O that treat the second and third kind of terms 
appearing in the right hand side of the previous equation. 

Fact 2 : terms of the form 

with I < j < m are controlled in by the data. 
Fact 3 : terms of the form 

with I < j < m and <T < fh — j , are controlled in by the data. 
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The first kind of terms in equation (|6ip is reduced again moving Vt past of Vq. 
We compute 

Again the second and third kind of terms in the previous equation are treated using 
Fact 2 and Fact 3. We keep going hkc this until we get a last term to be treated. 
This term has the form 

(62) VaiV'^'+'n). 

We must treat this term following another route. We write 

where Pj,^ is the horizontal projection. We make the expansion 

v.(v-(n^'^))= y: v.((v?n''')(v^^)). 

a+/3— m 

Let us consider the expression ^T^b^ ■ We compute 

By Fact {J\ we know each summand is controlled in iJ^^^ and therefore the full 
expression is. Note by this, that if Vq is applied to them, we get by writing 
Va = Pa" Va' — ToV T ■, that thc outcome is controlled in Hl'^^ . Let us consider 
now the expression We compute 

(63) = Y (V^^)(V^V,,A^). 

1+8=13 

Note that if we expand V^l/iV (using the quotient rule) we get using the data that 
this term is controlled at least in Hlf^^ . We consider next the expression Vf,'A^. 
From the identity 

we compute 

(64) V^Vf,' = V^rf^Vb'VxN - V^^^n^^V^A). 

We use the next Fact [Fact^ to treat the second term in the RHS of the previous 
equation. 

Fact 4 • Say 0<^ — 5<m— 1 and < S < ffi — 1. Then, terms of the form 

v^-*"(n,fv,(v^A)), 

are controlled in Hlf^^ by the data. 
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If we keep moving Vt past of Vh' in the RHS of equation (IM)) and applying at 
each time the Facf[^ we get that, except for the term 

occurring when (3 — ifi, all the rest are controlled in H^^^. Similarly if we apply 
Vq over the expression we get, after writing Vq. = P^" Vq/ — To^t, that 
when P < fh, the outcome is controlled in H^'^^ and when /? = m it can be written 
as the term 

plus a term controlled in H^'^^ . Putting all together we conclude that the expression 
(p^ is controlled in if we can prove that the expression 

is controlled in H^. To see that write it in the form 

and use the data. To finish the proposition it remains to prove Facts [7][^ that we 

do next. 

Proof of FactUi 

We prove it by induction. First note that 'VxTa — T'^Ilca is controlled in 
Hl+'^+^ and thus the Pact Q] holds when S = 1. A ssume we have shown the Pact [7] 
is valid until S — Sq, we will show it is also valid when (5 = (5o + 1. We compute 

a+l3=5o 

From the IH we know V^II is controlled in i/,*^^™ ^''^^ and by the assumption 
V^T" is controlled in The two estimates imply the Pact [3 □ 

Proof of Fact[B 
We compute 

a+/3+7— 771— J 

By hypothesis is controlled in and by Fact [7] any one of the other 

two factors in the previous equation is controlled at least in Hf~^^. The Pact gets 
proved from both estimates. □ 

Proof of Fact\M 

The proof of this fact follows by induction. First note that the Pact is valid 
when r = 0. This follows from Pact[7]and writing Vc = P^,^ Vc' — TcVt- Assume 
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the Fact\3iis valid for all (r, j) satisfying r<ro, ro<m-2,0<r<TO-j and 
1 ^ J ^ "T,. Then we show the Fact is valid when F = To + 1 as well. We compute 

Q+/3=r 

The factors V^II are controlled in a,nd therefore controlled in Hlf^'^ . 

It is enough to prove then that the factors V^Vc(V^Tfc) are controlled in by 
the data. We compute them in the form 

V^(Vc(V^j,Tfc)) = V^^-\V,{V'+^n) + r™(V^j,T")Rm„,,b„ - H /»V™(V^j,Tfc)). 

The third kind of term in the RHS of the previous equation is controlled in 
by the data. The second kind is controlled in by the same argument that 
Fact [H was proved. For the first kind of term, we move Vt past of Vc again. 
This generates two new terms which as we have shown for the two last terms in 
the RHS of the previous equation, are controlled in by the assumption and 
Fact 2. We keep operating like this until we get a last term 'V c'^^-'T^. Writing 
Vc = Pc'^ —Tc'Vt- and using Fact[l\we get that this last term is also controlled 
in Hf, thus finishing the proof. □ 

Proof ofFact\^ 

We prove this Fact by induction. First note that the Fact is valid when 
6 = 6 for 6 = 0,. .. ,171- 1. This follows directly by writing Vc P/'^c' - T^Vt 
and using the data. Assume now we have shown the Fact is valid for all {5, 5) 
with Q<5<5<ifi— 1 and 5 — 5 — L, where < L < m — 2. We will show the 
fact is valid also when 5 — 5 = L + 1. We write 

(65) v^+Hn,?Vc(v^7V)) = (v?n,?)(v^Vc(v^7V)), 

Q+/3=L + 1 

The factors V^^n are controlled in i/^^*"™ ")+'^^ -^^e need to show the factors 
V^(Vc(Vf,Ar)) are controlled in hI+'^ . We write 

v^(Vc(v^7V)) ^ v^-^VcVf+^Tv - n/VrfV?,7V). 

The second kind of term in the RHS of the previous equation is controlled in H^^^ 
by the assumption. For the first term we move again a Vt past of Vc and use 
the assumption. We keep moving V^s past of Vc until getting the last term 
VcV^'^'^A^. It follows from the data that this term is controlled in Hp'^ (note that 
/3 + 5 < TO — 1). Putting all together we get that the expression ([65]) is controlled 
in by the data. □ 

Proposition 20 Say f < r < r/_|_2- For satisfying < i < I and 1 < j < I 

andi+j < I+l, ||J(Wj)||^, AB(o,f)) is controlled by (the data) \\Wj^k\\ H'+^rg,^ 

/or fc = 0, . . . , j, \\'^^T''^\\m+^+''{B{o,r))' /or fc = 0, . . . ,j, f, r and \\k\\L2i^B{o,r))- 
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Proof: 

The proof is based analyzing the inductive formula for the current 

(66) J(Wj) = ^(v^^n)"! * • • • * (v™=n)"= * n' * vWfe 

(67) +^(vj'in)"i *---*(v^'=n)"= *n'* v«,(r*Rm*w^), 

where the indices satisfy nj(l + mj) + I + k = j , with k < j — 1 for summands 
appearing on the RHS of (|80p and nj {1 + frij) + k + 1 + q — j ~ 1 for summands 
of the form (|ST|) . We treat first the summands appearing on the RHS of ([50]) . 
We will show that the terms V™n and 11 are controlled in iJ^^^ and that the 
term VW^ is controlled in H^. The product that each summand of the form ([50)1 
represents is then controlled in H^. By Proposition [191 VJ'TI is controlled in H^^^ 
from ||Wfc||^m+i+i-fc , k — 0, . . . ,m and W^^^^ N\\ j^i+a+k k — 0, . . . ,m as long as 
m + i < I (and to + i + 2 > 2 which is satisfied trivially). We get this condition 
from the hypothesis : indeed we have 1 + to < j , i+ j < / + 1 so we have i + m < I. 
On the other hand we have a priori control on ||Wfc|j^i+j-fc for k = 0, . . . , j and 
II V^~'^A^||^i+2+fc k — 0,...,j from the data, which covers the condition on 
and N required before. Similarly we have control on ||n|| and we know 
i+j + l>i + 2. Let us consider now the factors VWfc where k < j — 1. Writing 
VaW/c = Vo'W/c + Ta Wfc+i we get that the expression is controlled in by 
the data. The proof that summands of the form ([ST]) are also controlled in is 
direct from what we have shown and the data, after expanding V^(T * Rm* W^,) 
using the product rule. □ 

We state below the global versions of Propositions (THl (THl HZl HH and [^D] and 
whose proof is straightforward. That will be useful in the proof of Lemma [21 

Proposition 21 Say Y, is a compact H^^^-Riemannian three-manifold, where I > 

0. Then, \\Ric\\^i+i, ||i^||^i+2, and r/+3 are controlled by (the data) ||Wo||^i^, 

ll'^illif;^; ll-^IU^,' '^'^'^ ^1+2, where A is a H^^^ -canonic harmonic atlas. 

Proposition 22 Say E is a. compact H^^^ -manifold, where / > 0. Then, (the 
data) ||Wo||^7^, l|Wi||^i^, ||_ft'||^2, ||A^||l2 and r/+2 control ||A^||^/+3, where A is 
-canonic harmonic atlas. 

Proposition 23 Say Y, is a compact H^^^-Riemannian three-manifold, where I > 

1. Then, for any {i,j) satisfying < j < I and I < i < I, W^jWh' controlled by 
the data ||W,+i||^^-i , ||W,||^^-., ||Wo||^^-i, l|Wi||^^-., ||J(W^)||^^-i, ||i^||i2, 
V andrj^i, where A is a H^^^-cononic harmonic atlas. 

Proposition 24 Say Y is a compact H^^^ -Riemannian three-manifold, I > 0. 
Then, for any {m,i) satisfying m > 0, i > and 2<i + m</ + 2, || VyII||^i^ 

is controlled by || Wfc|| „i+(,Ti-fc)-i , for k — 0, . . . , to, WV^Nll „{m-k-)+i+i , for k = 
"a "a 

0,. . . ,m, ||7^||i2, V and ri+2, where A is a -canonic harmonic atlas. 
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Proposition 25 Say T, is a compact H^'^^-Riemannian three-manifold, where I > 
0. Then, for any (i, j) satisfying^ <i<I,l<j<I andi+j < I+l, |1 J(Wj)||^^, 
is controlled by (the data) ||Wj_fe||^i+fc , for k = 0,...,j, ||V:^ ''A^||^i+fc+2, for 
k — 0, . . . , j , \\K\\]^2, V and rj^2, where A is a H^^^ -canonic harmonic atlas. 

In the next proposition we make a last step before proving Lemma [2] We 
will denote with an upper-index (fc) the fc-th Lie derivatives in the time direction 
dt = NT. 

Proposition 26 Say Yi is a compact H^^^ -Riemannian three- manifold, I > 0. 
Then, for any (m, i) satisfying 2<i + m<I + 2,0<i,0<m (the data) V, ri^2, 
||Wfc||^i-i+(™-fc) , fork = 0,...,m, control || V^7V||^i+i+(™-fc) , 

for fc = 0, . . . ,TO, and \\K'^'^''\\ „i+(m-k) for k — 0, . . . ,m — 1 (if m ^ 0), where A is 

"a 

-canonic harmonic atlas. 

Remark 5 Note that the hypothesis made on the Weyl fields Wfe in Proposition!^ 
are the same as the hypothesis made for the Weyl fields inside the Proposition 
[24l Also note that we may have used the norms || VyA^|| j^fi+i instead of the norms 
II A^'^*^) ll^i+i as we can see from the identity N" = N'VtN that one set of norms 
control the other. Finally note that the conclusions extracted on the lapse in 
Proposition [521 are exactly the hypothesis on the lapse inside Proposition [Ml 

Proof: 

The proof proceeds by induction. We treat first the case when m = and 
2 < I < / + 2. Note that by Proposition [7| HATHls is controlled by ||Wo||^i-i, and 
V. Note also that standard elliptic estimates on the elliptic system 

{K)ijk — tij^Bim, 

V^'ATy = 0. 

show that ||i^||^i^ is controlled by ||Wo||^i-i, r/+2 and V. This proves the required 
control on K . Now, use this estimate for K on equations ([50)1 and ([51]) (that define 
£'(Wi) and i?(Wi)) to show that ||Wi||^i-2 is controlled by ||Wo||^i-i, r7+2 and 
V. Now use Proposition [21] to show that ||-/V||^i+i is controlled by ||Wo||jLfi-i, 
ri+2 and V. Finally from Proposition [21] we get that ||ffic||^i-i is controlled and 
therefore by Theorem [4] ||g||^i+i is controlled too. 

We treat next the case when m = 1 and 1 < i < / + 1. The fact that K is 
controlled in H^_^^ follows from Proposition 12T1 This in turn implies the estimate 
on g' = -'2NK. It remains to show that N is controlled in H^J'^ and N' in H^_^^ by 
the data. The first estimate follows from Proposition]^ For the second estimate 
we need to differentiate the lapse equation. For convenience we write it in the form 

-AN + {Rg + k^)N = 1. 
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Note from this form of the lapse equation that the TV is a function of the metric g 
only. Differentiating with respect to time {t — k) we get 

(68) - AN-+\K\^N-^ A'N +{Rg-+2k)N. 
For the derivative of the Laplacian we have [7] 

A7 = .g^^V.Vb/ - g'^H^aDi^'g-cb + big''' 9de)) ■ 
and for the derivative of the scalar curvature we have [7] 

Rg'= -Atrgg'+SgSgg'- < Ric,g'> . 

We can write then 

(69) AW = -2iVif"''V„V6iV - g°''(VaiV)(-2(V=iV)if,b - fcVfciV), 
and 

(70) Rg-^ 2fc(-l + N\K\^) - 2kiy aSIbN)K''^ + 2NkRiCabK''\ 

A direct inspection of equations (|69p and ([70)1 shows that the RHS of equation 
is controlled in i?^ by the data. The elliptic estimates of Proposition [T] (I) applied 
to the equation ([55)1 would show that N' is controlled in -ff^^ if we can show 
that the i?3i"norm of A^' is controlled by the data. To get that estimate multiply 
equation ([55)1 by N' and integrate (denote the right hand of equation ()68p as F). 
We get 

/ \VN]'' + \K\''{Nfdvg^ [ FN-dvg<{[ F^dvg)H[ {Nfdvg)"^. 

It is apparent from this that ||iV'||^i is controlled by the data. 

Assume now that the proposition has been proved for all (m, i) with m < fh 
and 2<i + m</ + 2. We will show it is valid when (m, i) — {fh + 1, i) with 2 < 
i + m+1 < / + 2 as well. Observe that the estimates we want to prove for the metric 
g follows from those we want to prove for the lapse N and the second fundamental 
form K using g = —2NK. Observe too, that the hypothesis when (m, i) = {fh+1, i) 
contains the hypothesis when (m, i) ~ {fh, i + As a consequence of this, we have 
automatic control over \\V'^N\\ ^i+rn-k+'i, for fc = 0, . . . , tti and \\K^'''>\\^i+fH+i-k, for 
A: = 0, . . . , m— 1 from the data (for (m, i) — {fh+1, i)). It remains to prove therefore 
that the data (for {m,i) = {ffi + l,i)) controls also ||A:(™)||^i+i and |1 V^+^A^H^s+i. 
We show first the control over K and then we show it over A^. It follows from the 
RemarkjSjand PropositionjMjthat || V^n|| ^n+i+i-k , for fc = 0, . . . , m are controlled 

A 

by the data. If Uab is a T-nuU and symmetric, space-time tensor, the time derivative 
of U has the expression 

(71) Uab = NVrUab + N{Uacnb' + C/facn,^), 
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It is easy to see that If is T-nuU and symmetric. Choose Uab = Kab = Pq"' Pb ^ '^a'b' ■ 
Using recm'sively equation (|7T|) together with the fact that || V^n||j:^A+7+i-fc, for 

fc = 0, . . . ,m is controUed by the data, it is direct to see that lliir^'")!! i+i is con- 

A 

trolled by the data too. 

We prove now that || V™^^Af||^i+i is controlled by the data too. This will come 
from differentiating the lapse equation m + 1-times. We will rely on the following 
fact whose proof will not be included and is straightforward 

Fact 5 : the fh+l-th time derivative of the lapse equation has an expression of the 
form 

where F is an expression controlled in Hj^ by the data (for {m,i) ~ {in + l,i)). 

With respect to the proof of this fact, we only mention that it can be proved 
differentiating equation ((68)) m-times, using the expressions ([69)) and (|70)) for the 
RHS of (|68)) . The derivative of the Ricci tensor is given by [7] 

Ric = ^ALg' - S*g{6gg) - ^Dgd{trgg'), 

where is the Lichnerowicz Laplacian (see [7]). The time derivative of the con- 
nection V is calculated by the following formula. Denote V* the covariant derivative 
of g{t) and pick an arbitrary time independent vector field 11°'. Then 

(72) V*C/^ = V^•f/''+f;;^[/^ 
where T is 

(73) rl = liVl'g{t)bd + ^l'git)ad - K'9(t)ab)g{tf^. 

Now Va' at t = is calculated by differentiating (|75)) with respect to time and 
evaluating at t = t^ (note that the time derivative of the first term in (|72p vanishes). 

From FaciEjand Proposition[l]we get therefore that |1 V™^^iV||^i+i is controlled 
by the data. □ 

We are ready to prove Lemma [2] 

Proof (of Lemma [2)): 

We prove first that the 1 + 2 harmonic radius is controlled from below by the 
BR-functional \\{g, K)\\br- We consider the case 7 = first and then the cases 
/ > 0. 

Case I — 0. This case follows by Proposition [8] and Theorem |4| 
Case I > 0. The proof in this case proceeds as follows. We prove first that 
'"7+2 is controlled by rj^i and ||(.g, -ftr)!!^^,. Indeed, the proof of that shows, more 
generally, that for any 1 < J < /, rj+2 is controlled by rj+i and \\{g, K)\\br- 
As a consequence, we deduce that rj^j is controlled by r2 and ||((7, fC)!!^^;. By 
case I = 0, r2 is controlled by ||(5, iC)||Bfl. The Case I > would then follow. 
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Figure 1 : Iteration of control of the Sobolev norms of Wq and Wi . Each column 
contains the kind of norms that control the norms on the previous column. In the 
table 1 = 1-1. 



We prove now that is controlled by rj^i and \\{g, K)\\b]i. Note that by 

Proposition [51 H^Hl^ is controlled by the BR-functional. We can then replace 
by \\{g, K)\\br in the hypothesis of Propositions and that we are 
going to use in what follows. By Proposition [5T] we know that rj_^2 is controlled by 
||Wo||„/-i, ||Wi|| /-I, rj_|_i and ||((7,i^)||_B_R,, where .4 is a /f^"*"^ -canonic harmonic 

A_ A 

atlas. If /— 1 = we are done, as we have that is controlled by r2 and \\{g, K)\\br- 
If instead / > 1 we apply consecutively Propositions [53] and [551 From Proposition 
[23] applied with I = I,i = I— 1 and j = 0, 1 we deduce that ||Wo||^j-i and 
||Wi||^r_i are controlled by ||Wo||^/-., ||Wi|i^r_., ||W2||^/-2, ||J(Wi1||^f_., 
rjj^i and ||((7, i4r)||Bi?- From Proposition [25l applied with 1 = 1— 1, j = l and 
i = 1 — 2 we deduce that || J(Wi)||^/_2 is controlled by ||Wo||^/-i, ||Wi||^/_2 and 
II {g, K)\\br- To organize visually the iteration, we have included the Figure[Tl The 
norms included at a given column control the norms at the previous column (there 
may be more information at a given column than what is actually needed to control 
the norms on a previous column). The step we have just done represents the control 
of the second column over the first. Now, applying consecutively Propositions 1231 
and 1251 over each new column we get at each step, we reach a last column with the 
norms ||Wo||//o^, . . . , ||Wj||^o^0. As a result of this we get that rj_|_2 is controlled 
by J-j+i and \\{g,K)\\BR as desired. 

To show control over \\K\\ „/+i observe that as we have shown above ||Wo|| „/-i 

"a _ "a 

and ||Wi|| „/_i are controlled by \\{g, K)\\br. where is a _ff^~^-canonic harmonic 
"a 

atlas. By Proposition [5T1 ||i^||„/+i is controlled by \\{g, K)\\br- It follows that 

"a 

||-ftr||„/+i is controlled by \\{g, K)\\br too. □ 

^^Note that f'/^x' ^ ^^"^ control ||Wj||^o for j = 0, . . . , I and where ^ is a H^+^-canonic 
harmonic atlas 
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To finish this section let us prove a proposition on the structure of the L^-norm 
of the current J(Wi) that will be of use in the initial value formulation. 

Proposition 27 Let T, be a compact three-manifold, {g, K) a cosmological nor- 
malized state and j > 2. Then the L^-norm of the current J(Wj) can be bounded 
as 

(74) ||J(W,)|U2 < C(£,_i,j.,V)q| +C(£:,_i,i.,V). 

Proof: 

The proof is better divided in two cases, when j = 2 and when j > 3. Before 
going into the analysis of these cases let us make the following observation. 

Observation 1. As an outcome of the proof of Lemma [2] it can be seen that 
llVJ'niU. form<j-2, ||Vr'n|Ui, |iVWo|Ui, ||VWfe||HO for fc = 0, . . . , j - 2 
and II J(Wj_i)|| jjo are controlled by v and V. 

Case j = 2. When j = 2 we can write schematically (see for instance Proposition 

j(W2) =n * n * vwo + Vt * vwq + n * vwi + t + Rm * Wi 

(75) 

+ n*T*Rm*Wo. 

We need to prove that the i^-norm of any one of the terms on the RHS of the 
last equation can be bounded by an expression of the form ((74|) . Let us treat 
them case by case. By Observation 1 in the first term in the RHS of equation 
((75|) the coefficients 11 of VWq are controlled in and the term VWq itself is 
controlled in iJ" by £i, v, and V. Altogether the L^-norm can therefore be bounded 
by C{£i,iy, V). Let us treat next the term VtH * VWq. We can writJ^ 

IIVtH* VW0IU2 <c(£i,j/,v)||VTn||Hi||vWo||^^i. 

By the Observation 1 the norm || VtHH^i is controlled by £1, v and V. We need 
to estimate the norm ||VWo||^fi. Write 

VaWo - P/ Va'Wo - T, Wi. 

Using formulas pT|l - ((26|) we can write schematically 

Wo^€*€*E + €*B*T + E*T*T. 

Using this and the connection formulas we compute (schematically) 

P/Va'^o =K *e*e* E + e*e*WE + K *e* B*T + e*WB *T + e* B* K 
+ WE *T *T + E *T * K. 

^"Note that for any [/ an V wc have \\U *y |j£2 < ||C/||/,4 || V||^4 , and still by Sobolev embeddings 
less or equal than C{£i, u, V)||y|| ji^i ||J7||^i . 
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From this we can write 

(76) ||P/V,,Wo|U. < C{£^,v,V){\\Eo\\hi + ||i?o||i/^). 
In the same way one can prove 

(77) \\VW4hi<C{S,,v,V){\\E4hi + \\B\\hi)- 

From the eUiptic regularity of Proposition[l]apphed to the elhptic system ([28|) - (|3T1) 
we have 

(78) \\{E^,B,)\\Hi<C{E,,vymE2,B^)\\H0 + \\{E^,B,)\\H0 + \^^^ 
and 

(79) \\{E^,Bo)\\hi <C(£i,y,VmE^,Bi)\\Hi+\\{E„,B^)\\HoJ. 
It follows from equations ^T}, and that 

II VtH * VW0IU2 < C{£i,v, V)Q| + C{£u V). 

Let us treat consider now the term 11 * VWi in the RHS of equation (|75|) . The 
coefhcient 11 is controlled in Hg by £1, v and V. Now use ([75)) to get a bound on 
||n* VWi||^2 of the desired form. The L^-norm of the last two terms r*Rm*Wi 
and n * r * Rm * Wq are treated similarly. 

Case j > 3. This time we appeal to the full structure of the current J(Wj) 
provided by Proposition [151 Recall 

(80) j(Wj) = ^(v?,'^n)"i * • • • * (v™=n)"= * n' * vWfe 

(81) + ^(v*in)"i * • • • * (v*=n)"= * n' * v^(t * Km * w^), 

The sum on the RHS of ([50)1 is among sequences ((mi, ni), . . . , (ms, ^s), I, k) with 
k < j — 1, mi >!,..., TTis > 1 and '^j(l + TT^j) + ^ + fc = j, while the sum ()8ip 
is among sequences ((mi, ni), . . . , (mg, ris), q, fc) with mi > l,...,ms > 1 and 
J2j nj{l + rhj) + k + I + q ^ j — 1. 

Let us treat the terms on the RHS of equation ([80)1 . Say k — 0. In this case 
the coefficients in front of VWq are controlled in Hg except when there is only 
one coefficient of the form V^^^II. Now by Observation 1 either VWo or V^^^II 
are controlled in Hg by £j-i,v,V. It follows that the norm ||V:jr^n * VWoHl^ 
is controlled by fi, v and V. Say now that k = 1, . . . , j — 2. In these cases the 
coefficients in front of are controlled in by £1, ly and V. It follows from 

this and Observation 1 that when fc = 0, . . . , j — 2 the L^-norm of the expression in 
the RHS of equation ()80p is controlled by £j^i, v and V. Say finally that k = j — 1. 
In this case there is only one possibility, namely the expression on the RHS of 
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equation ([50)1 is of the form 11 * VWj_i. From elliptic regularity applied to the 
elHptic system p8)) -([3T |) we know 

< c(£i,j.,v)(||(i?,,i?,)|Uo+|i(ii;,_ 

It follows that 

II VW,_i|U2 < C7(£i, J., V)(g| +C(£j_i,z.,V)), 

as desired. This finishes the treatment of the terms of the form ([SO)) . The terms of 
the form ([8T|) are treated similarly. □ 

3.4 Controlling the flow (g, K) along evolution. 

In this section we introduce the functional space in which solutions to the Ein- 
stein flow equations will lie. In the same vein as Section [3^ we introduce a BR- 
functional (see later) 

k=i-2 

||</'*g||si? = \\9\\c°(ia)(HA) + \\K\\c»(hl)(HA) + X! \\^^k,Bk)\\cO(Ifi). 

k=0 

and discuss how it controls: i. the set of flow solutions, ii. a dynamical smooth 
structure on the space-time manifold, and iii. the space-time metric as a metric 
over the space-time manifold with the dynamical smooth structure. The results 
will be used in the next section when we discuss the initial value formulation for 
the CMC gauge. 

We start giving the definition of admissible gauge. We then give the definitions 
of space-time solutions and flow solutions. After that we introduce the BR-norm 
and present the main results of the section. 

Definition 1 Say I is an interval (for the time coordinate t — k). Define 
C{I,a,(i){Hi,) = n^~Q (/, -fff ~"') where the subindex -k indicates the structure 
with respect to which the Sobolev space is defined. Each space C^{I,H^ "') is pro- 
vided with its usual sup norm. We may allow (3 to be less than a. 

Definition 2 Say {Y,,Aoo) is a C°° -manifold. A differentiable function 
X : C°°(5(S)) X C°°(r5(S)) C°°(r(S)) which is diffeomorphism invariant i.e. 
X(4>*g, (jy* K) = (j)*{X{g, K)) for any diffeomorphism : E — > E (the classes C°°(*) 
are defined with respect to Aoo) is an admissible shifl iff 

1. for any H^^^ (J > 2J atlas A which is compatible with Aqc at the j + 1 level 
of regularity, X can be extended uniquely to a differentiable function (also 
denoted by X) X : iJ^ x iJ^T^ ~^ property that: 

2. \\X\\^j is controlled by \\g\\^i and \\K\\^j-i and 
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3. for every path {g,K){t), such that {g' , K') is in C{I,j — k,j){Hj[) x C{I,j — 
k,j — 1){H^) (0 < k < j) and with the norms \\{g, K)\\c"{i.j)xC"{i.j) '^i^-d- 

\\i9',K')\\c(ij-k,j)xCii,j-k,j-i) bounded by A, X' is in C{I,j - k,j){HA) 
with norm controlled by A. 

Examples. 1. The zero shift is an example of an admissible gauge. 

2. Definition [2] of admissible gauges share many of its properties with the 
Andersson-Moncrief gauge I5j on hyperbolic manifolds. Although we do not claim 
that the Andersson-Moncrief gauge is admissible we would like to explain how some 
of the characteristics of admissible gauges are actually present in the it. Let gH 
be a hyperbolic metric on S. For every g in T, perform a diffeomorphism (j) in T, 
in such a way that the identity is a harmonic map between (j)*{g) and gn- The 
identity is harmonic iff the vector field V'' = g'^^e^lViCi — Vf^e^) is zero. The 
CMCSH (Constant mean curvature - spatially harmonic) fixes the shift X in such 
a way that at every time the identity map id : {'E,g{t)) — ^ CE^gn) is harmonic. 
With this condition the equation for the lapse and shift are the following [5] 

(82) -AN+\K\^N^1, 

AX^ + ffic'^X-^ - V^' =(-2Ari^"" + 2V"X")e'(V™e„ - V^e„) 
^ ' + 2V™iVi\:\„ - V'Nk. 

Thus, when the identity is a harmonic map, X is defined through ([SS)) and ((83)) 
and otherwise it is defined by making it diffemorphism invariant. Consider the C'°° 
atlas Aoc for which gn is C°°. As the differential operators ([82]) and ([83| defining 
(TV, X) are elliptic its is clear by elliptic regularity that if {g, K) are C°° then {N, X) 
are too. Therefore X : C°°(5(E)) x C°°(T5(E)) C°°{TS{E)). In the same way 
if .4 is a atlas which is compatible with Aoo at the j-level of regularity, by 

elliptic regularity X extends to differentiable functions on 

the other hand Lemma 3.2 in [5] shows that is controlled by ||g||jLf3-i and 

||i4r|| j:^3-2 . This says that X satisfies item 2 in Definition [2] (in fact it represents an 
improvement). Also from Lemma 3.2 in [5] we have that, for a path (g, K){X), the 
norm |lX'||„i is controlled by ||5||rrj-i, ||i^||f,3-2, ||.g'|| rrj-i , \\K'\\„j-2. The shift 

A A I I A A ^A 

X satisfies therefore item 3 in Definition [21 with k ^ j (in fact it represents and 
improvement). 

Remark 6 Observe the following property of admissible gauges. Say (S,-4) is a 
C°°-three-manifold and suppose X is admissible. Suppose Ai is another C°° atlas 
in E compatible with A at least at the third level of regularity. It is known there is 
a C°°-diffeomorphism (j) : (S,-4) {I],Ai) and therefore the properties of X over 
A pull back over Ai . 

Definition 3 A W -space-time solution (M,g) of the Einstein equations (in vac- 
uum) is a H^^^-Lorentzian four-manifold (M,g) satisfying Ric = 0. 
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Definition 4 A H^-flow solution oj the Einstein CMC flow equations in the admis- 
sible spatial gauge X{g,K), is a space-time solution (M,g) (in vacuum) for which 
there is a H'' -diffeomorphism : M — > E x J where (E, A) is a H^^^ -three-manifold 
and T, X I is supplied with the product structure, such that 

1. 4>*{g) is (therefore) a H^^^ -space-time solution of the Einstein equations in 
vacuum. 

2. The time foliation is CMC. 

3. The components g, N, X of the 3 + 1 splitting of (j)^g in T, x I satisfy: g is in 
C{I,W_^^) andg' is m C{I ,i-l,i-l){HA). N,X are mC{I,i-l,i-l){HA), 
and K is in C{I,i — l,i — 1){Ha). These fields moreover satisfy the CMC 
flow equations (Q)-®- 

4-. For every k — 0, . . . i~2, (/)*(Wfe) has electric-magnetic decomposition {Ek, Bk) 
in C{I , z — 2 — fc,i — 2 — k){Hj(), satisfying the equations \2^) - ^3l\) . 

5. There is a dynamical H^^^ -atlas A{t) with dynamical charts {x^{ty\ for which 
the set of space-time charts {x^ = form a H^^^-atlas of T, x I 

making (j) a H^'^^ -diffeomorphism. Also the transition functions Xa(xp) are 
in C{I,i + \,i + \)(Hxi^) and {g,N,X:^^) are in C{I,i,i){Hx^), where X:^:^ is 
the .shift vector of the coordinate system . 

A comment is in order. In item 5 of Definition [H tlie time derivatives required 
to define the spaces C{I,-k,-k){Hxij) are witli respect to the time coordinate of x^. 

In practice (for i — 3) the dynamical atlas A{t) will be given out of the following 
construction of harmonic coordinates. Let {x''} be a coordinate chart in a H'^- 
Riemannian manifold (E, g). Say B{o, ai) C B{o, 02) are two balls inside the chart. 
Pick a smooth non-negative function ^ being one in B(o, ai) and zero in B{o, 02) 
and define for any given metric g in E, the metric gij — S,gij + (1 Think the 

coordinates {x'^} as coordinates on a three-torus T^ with metric g over the chart 
and extended to be flat on the rest. Extend the coordinate x'^ smoothly to the rest 
of the torus in such a way that Jj.^ x^dv^ where g is cjij = S^gij + (1 — O^ij- Define 
the function h^ to be zero on the chart and equal to l^gX^ on the rest of the torus. 
Then we can solve uniquely for x^ in ^gX^ = h^ if we impose the condition that 

x^dvg = 0. Now suppose A = {{xa}, a = 1, . . . , n} is a TJ^-canonic harmonic 
atlas for the manifold (E,g) where each {xa} is defined over a ball B{oa,r2/2) 
{r2{g)). Pick S < 1 but close to one. Choose ai = Sr2 and a2 — r2. Then for 
e sufficiently small, if \\g — 5||^2 < e the functions {x'^,k = 1,2,3} defined over 
B{oa, '"2/2) are harmonic coordinates for g. Also if e is sufficiently small, the charts 
{xa} extend to B{oa,Sr2) and satisfying 

i53 4 
(84) —Sjk < gjk < -^Sjk, 
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(85) J2 I \^9,k?dv.)<- 

Thus we have constructed a new harmonic atlas for the metric g that we will call S- 
canonic. It is clear that all the results we have discussed so far in Section [3. 31 hold if 
instead of using a canonic atlas we use (5-canonic atlas for a fixed i^ . Now suppose 
g{t) is a path of metrics in C°(/, 2)(i?_^) and suppose \\g{0) — < e/2 (with e 

as before), then for any t in a subinterval /' C / the harmonic atlas {xa{t)} = A{t) 
is well defined. We will call A{t) the dynamical atlas constructed out of (g. A, g{t)) 
and denote the dynamical charts by {xa{t)}- 

Remark 7 As constructed the atlas A{t) is only an iJ "^-harmonic dynamical atlas 
and not, as is required, a 7?'*-harmonic dynamical atlas. We will see later (Propo- 
sition [30|) that in fact this construction gives us the desired i?'*-harmonic atlas. 

We describe now a useful i3i?-norm (Bel- Robinson norm) on the space of metrics 
on E X / with the properties in Definition [4] above and then move to explain how 
it controls the space of flow solutions. 

Definition 5 Let be a flow metric on x I as explained in Definition 
Define the BR-norm 0/ g out of A by 

\\(l)*g\\BB, = \\9\\c°{I,2)(Ha) + \\K\\c'>{I,1)(Ha)+ X! II(^J'-^j)IIcO(/,0)(Ha) 

Note that the definition of BR-norm does not need the condition in item 4 of 
Definition m to be imposed on 0*g. The importance of the i?i?-norm is that it is 
easy to handle in the Einstein equations and is intended to measure (some) ^^H\xr 
norm" of 0*(g) without being in H\^j. At the same time as explained later it does 
controls (some) iJ^^^ -norms with respect to the coordinates Xq,. 

From now on we fix z = 3. We note however that the treatment we make may 
well be systematized to any regularity i > 4. 

We will need elliptic estimates for scalar equations of the form 

A(/) = h, 

where h has an expression of the form 

/i = To o Vb -I- Ti o VFi +T20 VV2, 

for tensors and Vm, rn = 0, 1,2 of arbitrary rank and where the operation o 
is any full contraction (like T^'"^VaVcV2,fc). In particular we would like to obtain 
elliptic estimates in terms of the L^-norm of Vi, i = 1,2,3. (see the hypothesis 
inside the Proposition). Thus the non- homogeneous term h will lie, in general, in 

^^We have said what a i?^-<5-canonic is. It should be clear how to define H*-<5-canonic. 
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a negative Sobolev space. In this context, the eUiptic estimates that we wiU obtain 
(except item 1 in the proposition below) do not follow from standard elliptic esti- 
mates or Proposition [Hand deserve a special treatmenj^ due to the low regularity 
of the metric g which is assumed in . 

Proposition 28 Let {Yi,A,g) he a C°° -Riemannian manifold and say A is a 
canonic harmonic atlas. Let T^, Vm, m = 0, 1, 2 he C°° -tensors of arbitrary rank. 
Say <f> is a solution of 

(86) A0 = h, 

with (j)dvg — 0. Then, we have the elliptic estimates 

1. For j = 0,1,2 we have 

U\\H^+2<C{r^,Vol)\\h\\^,^, 

2. Lfh = Tio VV^i + To o Vo, we have 

ll^ll^i < C(r3, Vol, llTill^^, llTollHi, WV^Wh'x, ll^olU^), 

3. Lfh = TooVQ+Tio WVi +T20 V2 we have 

Mh'X < C{r3,Vol, \\T,\\h^^, IITsIIh^, UToHh^, ^2!!//^, ^oIIh^). 

With low regularity the estimates extend in the following sense. Say {'E,,Ai,gi) 
is a sequence converging in LL^ (for Ai) and (for gi) to (I],^oo,<7oo) where 
(S, ,4oo, <7oo) is H'^-Riemannian manifold, with Aoo 0, -canonic harmonic atlas. 
Say too hi, Tm,i and Vm.i are sequences with hi — s- hoc in for j either 0, 1 or 2, 
Trn.i 7m. 00 in H'"^^^ (for item 2) and iJ™ (for item 3) and Vm,i —> Vm,oc in ■ 
Finally say (pi is the sequence of solutions to US 0(1 . Then, (for item 1) (f>i (poo in 
for j either 0, 1 or 2 (for item 2) (pi — > (poo in and (for item 3) pi ^ (poo 
in H°. 



Proof: 

Recall the fact that A : _ff ° is Fredholm of kernel the constants and index 

zero. Moreover if we consider the operators A : Hj^ iJ^ and A^^ : iJ^ 
as operators from the orthogonal complement of the constants into the orthogonal 
complement of the constants, their norms are controlled by r^ and l^o@. 

Item 1. This case follows from standard elliptic estimates (or Proposition [T]) 
and the observation above. 

^^We haven't found these estimates in the standard references on PDE. 

^''The boundedness of A is direct form elHptic estimates. The boundedness of follows 
by first proving the image of A is the orthogonal complement of the constants and then using 
the open mapping theorem. The fact that the image of A is the orthogonal complement of the 
constants follows by expanding in terms of eigenfunctions of A. 
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Item 3. Pick ^ with ^dvg — 0. Observe that in any contraction of the form 
T o WV or T o W^V we can always integrate by parts once or twice respectively. 
Multiply equation ([86|l by ^ and integrate. We get 

/ ^Acpdvg = / (v(eri)) o V, + {w^i^m o V2 + (eTo) o Vodvg. 

Expanding the covariant derivatives on the RHS of the last equation and using 
Sobolev embeddings and Holder inequalities we get 

^^^^ l^iAOHvg <C(r-3,^oO(||eilifi(l|Vo||H^||ro||ifO +||T2||^^||T/2||if^) 

+ U\\h^JT2\\hi\\V2\\h'X + U\\Hl\\Tl\\Hi,\\Vl\\H'x)- 

Make f = A^. Then from ||^||^2^ < C{r3,Vol)U\\H0^ we get that f ->< > 
is a bounded linear map from Lg into K. In particular ||(/)||i2 and so ||(/'||^fO are 
controlled by r^, Vol and the respective norms of Tm and Vm for m = 0, 1, 2. 
/iem 2. Multiply equation ([86]) by and integrate by parts. We get 

/ |V(/)pdwg = / \7(j)oTioVi+(l)S/TioVi+<t>noVo. 

Using again Sobolev embeddings and Holder inequalities we have 

(88) lj^(f>\'dvg<C{r,,Vol){MH^jT,\\H2jV^\\H<x + M 

From item 3 we have control on ||(/'||/fO^. Using it and equation (j88p we get control 
on 1 1 01 1 //I from r^, Vol and the respective norms of Tm and Vm for to = 0, 1, 2 as 
desired. 

We prove now (sketchily) the last part of the proposition. Assume then we have 
a sequence 5^, 0^, /i^, Vm,i, Tm,i as indicated in the statement of the proposition. 

Item 1 Although the proof of this case is trivial, we will follow one that can be 
applied in the proofs of items 2 and 3 too. For a metric g and a field write 

= Ag^^ + g'^^Tl.V,^ + ig-' - gt)VaVb<P. 

where V - V°° = F. Note that ||F,||^2^ ^ and \\gf - g^Wn^^ ^ 0. Using 
this formula rewrite the subtraction of equation (ISS)) for the fields <t>i,gu hi to the 
equation ((86)) for the fields (j)k,gk,hk, as 

with ||/i;A;||/fj and (pi — (jy^dyg^ — > 0. It follows that — (/)fc||^i+2 — > 0, 

thus the sequence {4>i} is Cauchy in H^^^ and therefore convergent. 

Items 2 and 3. This item follows in the same way as item 1 was proved. Observe 
that the RHS of equations (|87p and ([55)1 tends to zero (thus making the norms 
11011^0 and tend to zero too) if, simultaneously, the (respective) norm of 
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one of the fields in each one of the pairs (Tq, Vq), (Ti, Vi) and (T2, V2) tends to zero. 
In the same way as in item 1 subtract equation (j86p for the fields 4>i, gi, Tm,u ^711,1 
for m = 0, 1, 2 to the equation ([86| for the fields (f)k,gk^ Tm,k, Vm,k for m = 0, 1, 2. 
Rearranging the result conveniently and using the observation above we get that 
4>i — (/ifc is a Cauchy sequence in (for item 2) and in (for item 3). The 

claim then follows. □ 

Remark 8 One can obtain exactly the same kind of estimates as in items 1 and 
2 in Proposition for equations of the form 

where / is in and / > > 0. This time the estimates depend also on the 
ijji-norm of / and fmf- The estimates, in particular, will be required to estimate 
the i?^-norm of A^'" (the third time derivative of the lapse). 

Proposition 29 Let be a flow solution in x I with BR-norm out of a H'^- 
atlas A hounded by K. Then, the norms of g, g' and K in C^{I' , H'^), C(I' , 2, 2)(iJ^) 
and C{I' , 2, 2){Hj\) respectively and the norms of N , X and X' in C{I' , 3, 3){H^), 
C'^{r ,2){Hj[) and C(/', 2, 2)(i?^) respectively, are controlled by A, where I' is a 
subinterval of I with size controlled from below by A. 

Proof: 

The idea for the proof is to look at the equations 

(89) g-^-2NK + Cxg, 

(90) K-=-W'^N + N{E-KoK) + CxK, 

(91) -AN+\K\^N = 1, 

(92) E- = NCurlB - {\/N A B) - ^N{E x K) ~ ^N{E.K)g - hzNE, +CxE, 

(93) B' = -NCurlE ~ {V N A E) - ^N{B x K) - ^N{B.K)g - ^kNB + CxB, 

and observe that the time derivatives of g and K can be calculated using recursively 
those equations and its time derivatives. There are delicate points however that 
have to de addressed with care. 

Without loss of generality we will prove the proposition with A a canonical atlas 
for the initial metric g(oj^. We will make that assumption on A from now on. 

Observation 1. By the definition of the _Bi?-norm, g is in C°{I,2){H_a) with 
norm controlled by A. Therefore, given e > 0, there is an interval I' C I with a 

^■^The initial metric g{0) is always assumed to be in H^^ 
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size controlled from below by e and A such that the atlas A{t) is well defined for 
all t £ I' and every chart satisfies ||a;J;(i) — x^{0)\\c{i' ,3){h^) < As a result, if e is 
chosen sufficiently small, then for any t £ I' the identity id : (E,^(0)) {^,A{t)) 
is a ff'^-diffeomorphisrro. Moreover, for any tensor field U, the norm llC^II-ff^^jj 
controls the norm ||C/||_ff^ and vice versa, for any t e I' . 

Observation 2. As an application of the Observation 1, note that the norms 

llffllc«(/,2)(ff^(t)), \\K\\c"{I^)(H^^t)), \\{Eo,Bo)\\cO{IM){H^(f,), ll(£l^l)llcO(/,0)(H^(t)) 

are controlled by A. Therefore, by the results in Section 13.3.21 for any t E I' 

the norms bllffj^^^,, H-^H^iw' H^l'^iit,' IK^o, ^o)!!//^^^, , || (-Bi, Bi)||//o^^^^ are con- 
trolled by A. From the Observation 1 above so are controlled the same norms but 
with respect to the atlas A. 

The proof goes in a ladder-like scheme. We follow the order 

g,K,N,X^g\K\N\X-^g-,K-\N-\X-^g-;N-\X-. 

Namely, proceeding from left to right, once we prove that a field is controlled (in 
its respective norm) we use that control to prove control on the next field (in its 
respective norm). 

1- \\9\\c''(r .2)(Hji,) controlled by A. Direct from the definition of BR-norm. 

2. Il-ft^ ||c''(/'.2)(_ff^) controlled by A. Direct from Observation 2. 

3. ||-/V||co(/',3)(ff^) controlled by A. Direct from Observation 2. 

4. II Ar||(70(/'^2)(ff^) controlled by A. As X is admissible we get, applying item 1 in 
the Definition [2] with j = 3 and with atlas A{t)^ that ||X||^3^^ ^ is controlled by A. 
This item (4) then follows by Observation 1. 

5- Wg'Wcii' ,2){Ha) controlled by A. From equation (|89p and Observation 2 (or items 
2 and 5) we see that we need to control the iJjj-norm of Cxg- This follows from 
the formula 

{Cxg)ab = '^aXb + VbXa, 

and Observation 1. 

6. ||-ft"||co(/',i)(H^) controlled by A. From equation (|90p and Observation 2 we see 
that we need to control the Hj^-noim of CxK . This follows from the formula 

(94) CxUab = VxUab + Uac^bX' + UbcVaX^. 

with U = K, Observation 1 and item 4- 

7. ||-/Vl|c°(7',2)(//^) controlled by A. Differentiating the lapse equation with respect 
to time we get 

(95) - AN'+\K\'^N- = -{\K\^yN + AW. 
Recall that the derivative of the Laplacian has the expression 

(96) A7 = s^'VaVb/ - g''\Vaf)iV'g-cb + l^big'^gde))- 

^^The fact that id is a difTeomorphism is somehow standard and can be elaborated from the 
following fact. Say (/> : — > is a H^^ map being the identity on -8(0, 1)'^ and with ||0 — 
^'^llfrs 5: ^- Then if t is sufficiently small </>~^ exists, is in H^^ and has _ff^3 norm controlled by e 
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Using (IM)) (with / = N) in ([TO)) and the previous items, we see that the RHS 
of ([95]) is controUed in -ff3l(t) • '^^'^ elliptic estimates of Proposition [28] (see the 
Remark[8|) show that N' is controlled in (in fact in ^^^(t)) by A and therefore 

by Observation 1 so in H^. 

8. ||^'||c''(/',2)(_ff^) controlled by A. As X is admissible we can apply item 2 in 
Definition [2] with j = 2, k = 2 and atlas to get that X' is controlled in by A. 

9- II5 "llco(/'.i)(-ffA) controlled by A. Differentiating equation (|89p in time we get 

(97) (7--=-2Ar-i^-2Ar(i^)-+(£x.9)'. 

The first two terms in the previous equation are controlled in by the items 6 
and 7. Let us consider the derivative in time of the Lie derivative of g with respect 
to X. Differentiating Cxgab — ^ aXb + ^bXa with respect to time we get 

(98) {Cx9)ab - ^a'Xb + VfeX, + WaXl + VbX\. 

Recalling how the time derivative of V was calculated in Proposition [26l a direct 
inspection of all the terms of VaXf, show that Va'Xb is controlled in H\. The 
control in H\ of the last two terms in equation ([98]) follows from item 8. 

10. Il^i^ "||c''(/',o)(i/^) controlled by A. Differentiate equation ([50)1 with respect to 
time. We get 

(99) K"^ ~-{V^yN - V^N'+NXE - K o K) + N{E'- {K o K)) + (CxK): 

A direct inspection of this equation and use of the previous items shows that all 
the terms in it, except perhaps NE' and {CxK)\ are controlled in by A. To 
check that the term NE' is controlled in examine equation Q and use the 
previous items. Recalling that HxK = VxKat + Kac^b^'^ + KtcVaX'^ and using 
the previous items, we see that the term {CxK)' is controlled in by A. 

11. II A^ "||co(/'.i)(_ff^) controlled by A. Differentiate (j95p once more with respect to 
time to get 

(100) - AN"+ \K\'^N"^2A'N'^ 2{\K\'^yN' + A"N + {\K\'^)"N. 

We want to apply Proposition [28] (see Remark [8]) . We check that the RHS of 
equation (jlOOp is controlled in by A. The first term AW in the RHS of 
equation (|100[) is calculated using equation ([96]) with f = N'. A direct inspection 
of its terms using items 5 and 7 show that they are controlled in by A. Now 
let us consider the time derivative of \K\'^. Writing lii'p = K abK ^d.g'^'^ g^'^ , the time 
derivative is 

(101) \K\'''^2KabK,dg-^-g'"' + 2K'abK,dg'"'g'"'- 

Differentiating once more we get 

\K\^"^2K"abK,d9^'g"' + 2K'abK',d9-^-g'"' + iK'abK.dg-^'g'"' 
iK'abK.dg'^'g"' + 2KabK,dg-^'g'"' + 2KabK,ag-^^g'"'. 
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Inspecting equations (|10ip and (|102p and using items 5,6,9 and 10 one gets that 
\K\'^' and are controlled in H^^ and i?^ respectively by A. This shows that 

the second and fourth terms in equation (jlOOp are controlled in Hy^. Finally, let 
us consider the term A"N. For this we differentiate A' in equation ([96|) once more 
(leaving f — N invariant). We get 

+ g^'VaNiV-'g-cb + V^g-'cb + ^^bi^g-^'g-de + g^'g-'de))- 

A direct inspection of this equation using items 1, 3, 5 and 9 and the expression of 
V found before shows that A"iV is also controlled in Hj^ by A. 

12. 11^ "||c"'(/',i)(ffyi) controlled by A. As X is admissible this item, follows applying 
item two in Definition [2] with j ^ 2, k ~ 1 and atlas A. 

13. To control ' l|c«(/',o)(ff^) , '\\c°{r ,o){,Ha) and ||X"1|^q^ proceed in the 
same fashion as in the previous items. □ 

Let {xa} be the (5— canonic dynamical as was defined after DefinitionHl Note as 
it was defined, it is only a i/'^-atlas and not a iJ'^-atlas as required in Definition 01 
In the proposition below we discuss this and further properties of A{t) in relation 
with Definition 21 Recall the coordinates defined over B{oa,r2/2) where 

is the i?^-harmonic radius of the metric g. In addition those charts can be extended 
to B[oa,5r2) satisfying the properties ([84]) and (fSS]) . 

Proposition 30 A flow solution (/)*g in x I with BR-norm hounded by A has 
the following properties (at least on a smaller interval I{A)) 

1. A controls the C{I ,A,'i){Hxfj) norm oj the transition Junctions Xa{xfi). 

2. A controls the C{1 ,3,3){Hxij) norm ofg~ (g, A^, ATx^) where ATx^ is the shifl 
vector of the coordinates X/3 . 

3. Let g(A) be a path of metrics parameterized by A. Then the BR-norm of 
g' = If 

j=i 

ll(</'*g)'ll = ll5llcfH/.2)(if^(o)) + ll^1lcfH/a)(Jf^(o))+Zlll(^i'^j)llc''(/,0)(ff^(o))' 

controls the norms ofx'^ — dxa/dX and g' in C{I, 4, 4)(_ffj;^) and C{I, 3, 3){Hxi^^) 
respectively. We will assume without lost of generality that the BR-norm of 
(0*g)' is bounded by A too (i.e. the same as the bound for the BR-norm of 
the solution). 

Remark 9 i. In item 1 the time derivatives oixa (needed in the norm C{1 , 4, 4)(_ffj.^)) 
are with respect to the time vector field defined by the chart {a;^(t),i}. 
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ii. Note that item 1 in Proposition 1301 implies that when E x / is provided with 
the set of charts t)}, a = 1, . . . , to} it makes it a i/''-fom'-manifold. 

iii. Item 2 in Proposition [30] shows that g is a /f^-metric in S x / with the 
{{{xa{t)Tt)}} atlas. 

iv. In item 3 we use the convention that (as we have been assuming) the 
C(/, 3, 3)(iJa;3)-norm of g' is the C(/, 3, 3)(7Ja;^)-norm of each of its horizontal 
components. To find which the components are, pick two A-independent horizontal 
vectors and and compute 

and 

Observing that 

d 

= ( — )' = N'T + NT' + X', 
^dt' ' 

we get 
Finally 

N' 

g'^.T^T" = -2gafcT"T"' = 2—. 

It follows that the components of g' form a set equivalent to {g' , N' , X'). If we 
prove that {g' , N' , X') arc in C{I,3,3){Hxij) with norm controlled by A the same 
will be true for g'. We will do that when proving the proposition. 

V. It is straightforward to prove from item 3 in Proposition [501 that {g',N',x') 
is in C(/,3,3)(i?^) and {g';N';X') is in C{I,2,2){Ha) with norm controlled by 
A. We will use this fact when discussing the initial value formulation in the CMC 
gauge. 

Proof: 

Item 1. Step 1. We show first that the transition functions Xo,{xp) have 
H^^ norm controlled by A. By Proposition [21] we know \\Ric{t)\\ui^ is controlled 
by A. For the same argument as in Proposition [2] this implies that the norm 
llsll^fl {B(oi},Sr2)) controlled by A. Now, express the harmonic condition of the 
coordinate xjj over {xp}. We have 




We can apply the elliptic estimates of Proposition [T] (I) to get that 
\\xa\\H^^{B(o^,,r2/2)nB{o^,r2/2)) IS Controlled by A. 

Step 2. Next we prove that the norms \\xa\\Hl^(B(o^,r2/2))-, \\xa"\\Hl^(B{o^,r2/2))-, 
\\^°'"'\\HUB(o^,r2/2)) and \\xa""\\H»jB(o^,r2/2)) are controlled by A and where the 
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dot means derivative with respect to dt — NT + X (and not the time vector field 
of the chart {(x/3(t),i)} as we will do in Step 3). We will prove that studying the 
defining equation for on the torus. Control will be obtained in the order 

We have defined by solving 

(103) A,(,)x^ = h", 
subject to the condition 

(104) / xldvg^,)^0. 

The idea is to differentiate consecutively these equations with respect to time an in 
each step use the control gotten in the previous step. Differentiating both equations 
with respect to time we get 

(105) A(x^)=-Axt 
and 

(106) / {xiYdv-g^t)^ - f xia-Nk + W'Xajdvg. 

Using formula (|96p for the derivative of the Laplacian we see that the RHS of 
equation (|105[) is controlled in H^^ by A. We split (a;Jj)' = (a;^)'o + c where c is 
a constant and {x'^)'o has average zero. From Proposition [28] we conclude that 
ll(2;Q)'o||_f/3 is controlled by A and using equation (|106p we conclude that c is 
controlled by A too. As a result ||a;Q'||/f3 (B(oa,r2/2)) is controlled by A too. 

We show next that the norms \\xa"\\H^^{B{o^,r2/2)), ||a;a"1|//i^(B(o„, 7-2/2)) and 
\\xa Who (B(o„,r2/2)) are controlled by A. We proceed exactly as we did before. 
Differentiate equations ()105p and ()106|) and use Proposition [28l A lengthy but 
straightforward calculation shows that every time equation (|105p is differentiated 
we can apply Proposition 1281 In other words the RHS of the equation 



71 — t 

Ax'<''> = ^ c„A(")x^('-") 



((★) denotes the ^-derivative) for each i = 1, 2, 3, 4 has the structure of h in Propo- 
sition The claim the follows. 

Step 3. We are now in condition to finish the proof of item 1. We will use the 
notation H*^ instead of H*^{B{op,r2/2) n B{oaipha,r2/2)). Let us denote xp = 
(xp(tf}),tf}) and Xo = (x/3(0),to) where {xq = a;^(0)} is a time independent chart. 
The functions tfj and t^ are the same and just t (the subindex is to differentiate dt^ 
from 9t^)- For any function / the chain rule gives 

df df , df dxl 



dto dtp dxp dto 
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From this we conclude that 

d dxi 
(107) —=NT + X- ^ 



dtp dto dxy 

Now suppose we want to show that H^taa^Qll^fa is controlled by A. Apply dtp using 
formula ^07]) over (for fc = 1, 2, 3). We get 

_ k- _ i ■ 

dtp " " " dxy 



Steps 1 and 2 show that both terms on the RHS of the previous equation are 
controlled in . Showing control over the other time derivatives follows the same 
strategy. To illustrate the procedure let us prove that ||9j^„a;J^||i/2 is controlled by 
A too. The remaining two cases are done in the same fashion. Apply the expression 
P07p twice over x^. We get 

^^ a^a^l dx^_d_ dx^dxi 

dtl dtl dto dx'p ^ dto ^ dto ^ dto dx\^ ' ^ dto dx'^ ^ dto dx^^ 

From Steps 1 and 2 we see that all the terms in the RHS of the previous equation 
except perhaps the third are controlled in iJ^^ . To treat the third we compute 

d dx^a d^x^fidx'' Sii d^x'' 

^ ' dto^ dto dx^p ' " dtl dx\^ dto dtodx^p ' 

The first term on the RHS of this last equation is controlled in H"^^ by Steps 1 and 
2. Let us consider the factor 

(109) 

dto dxp 

Qj.^ 

on the second term of the RHS of equation l|108p . As we know that the factor 

is controlled in H^^ we would like to prove that the term (|109p is also controlled 

in Hl^. Rewrite ((lB9t as 

d dx^ dx'o . _ 9^2;™ dx™ dx^ d ,dx^ 
dh^dx^W~' ~ dxfdU)'dx^ ^ dx^dt^^'dx^'' 

The first term in the previous equation is controlled by A in H^^ by Steps 1 and 
2. The second term is treated using 

dx'o dx'p _ 



From this we compute 



dx'pdxi 



d'^x^ _ dx^dxl^ d'^x'p 



dtodx'-p dx'-p dx^f^ dtodx^ ' 
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It follows from Steps 1 and 2 that this terms is controlled in H"^^ too. 

Item 2. Checking this property is straightforward. We will use the notation H* 
instead of-ff*^(B(o^,r2/2)). We will only check that ||Xx^||_f/3 and jj/Cat^Xx^ ||_f/2 
are controlled by A. The remaining two cases are carried out in the same fashion. 
From equation ll07l we get = X— gfr- We know ||X||/f3 is controlled by A 
and from item 1 we know ||9tgXo||/f3 is controlled by A. It follows that ||Xxg|j//3 is 
controlled by A too. Let us consider now the Lie derivative of Xx^ in the direction 

of dtp. We apply Cata first over ^j^^fr and then over X. We compute 

dto dx'f^ ^dtpdto'dxy 



and 



dt^'dt^ ~ ~dq " dta dxy^'' 



It follows from the equations above and item 1 that the 77^^ -norm of ^dt/i^^-^rr 
is controlled by A. Consider now Cdt^X. We compute 

Cdt X ~ Cq, X — C a^i X, 
"*I3 "'0 _A_o 

9*0 ax 



and 



dxl d 
£a^i X = —^\/ X + X''\/J^ 



^13 



dto ax^ dto dx^p ' 



It follows from these two equations, item 1 and Proposition that XWm 
is controlled by A. 

Item 3. We will show that x'p and g',N', X'^^ are controlled in H^^ and H^^ 
respectively. Here, prime means derivation with respect to A. The control of their 
time derivatives (as ^ or N'''') in their respective space, is carried out along 
similar lines and will not be included. We start showing that ||xo||//3 is controlled 

P x^ 

by A. We look at the equation defining {xp} in the torus. After differentiating in 
A we get 

(110) A(x^)' = -A'(x^), 
and 

(111) / (4)'d^;, = -/ x'p^^-^dv-g. 

As g' is controlled in Il\ by A so is controlled the term on the right of equation 
(jllip . Using formula it is seen that the RHS of equation (jllOp is controlled 
in iJy3 by A. The elliptic estimates of Proposition [^5] show that (x^)' is controlled 
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in H^^ by A. Still we want to show that |1 (a;^)'||/f4 is controlled by A. To do this 
step further we prove next that || (^iiO'Hifs is controlled by A. As a result if we 
write equation (|110p over the {xp} coordinates 

we can apply the elliptic estimates of Proposition [1] (I) to see that || (a;^)'|lif<i is 
controlled by A too. To show that H^'lj/fj is controlled by A we proceed as follows. 
We look at the system of equations 



(113) T^g'"'^l^„.g^J + Q{d.g*,9*\j = So,y + K^Ksj - kK,,, 



(114) V"£;o,ah = (i^ABo)fc, 

(115) {Curl Eo)ab = Bi^ab + \{Ba x K)ab - \kBo,ab- 

The first equation is on the {a;^(A)} coordinate system (the other two are tensorial) 
and for this reason we have to be specially careful when we differentiate with respect 
to A. We utilize the following. For any two tensor T, it is 

r.... r ^T(^ ^ \=T( ^ ^ 9x'"(0)fe'(0) 

^ ' ^ dx^ (A) ' dxi {\y ^ ax™ (o) ' dx^ (o) ' dx^ ( a) dxi ( a) ' 

Thus, when we differentiate with respect to A we get 

(117) (T.,)' = (n. +T.(^)' + r„,(^)'. 

Observe that (^^T(^)' is controlled in i?^^- Differentiating the system of equations 
pi3p - (lll5p with respect to A and evaluate it in the {xp\ coordinates. After a 
straightforward calculation we get an elliptic system of the form 

(118) g'^did^{g,,)' = {{E^%^Bl, 

(119) V\E^%^Bl 

(120) (CuTlE^)\,,,=Hl,, 

where H], is a term controlled in by A and iJ? and ffP are terms controlled 
in H^^ by A. Standard elliptic estimates show that ||(5ij)'||ff| is controlled by A. 

Now that we know x'^ is controlled in H^^ we conclude after inspecting equation 
pi7p that g' is controlled in by A. It is straightforward to deduce, after 
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differentiating equation ([5^ in A, that |lfC'||j:/2 is controlled by A. Differentiating 
the lapse equation and using standard elliptic estimates we get that N' is controlled 
in H^^. Using the definition of admissible gauge we get that X' is controlled in 
H^^ too. □ 

Before going into the inital value formulation let us mention that the space 
of flow solutions on an interval / is complete under the BR-norm. The proof is 
straightforward and will not be included. 

Proposition 31 Let (S,.4) be a H'^ -manifold. The space of H^- flow solutions 
over an interval I is complete under the BR-norm. 



4 Applications. 

4.1 The inital value formulation in the CMC gauge. 

Theorem 5 Let (S,^) be a H'^-manifold and say (go,i^o) 'is an initial state in 
X with ko <Q. Then 

1. There is a unique -flow solution over an interval L — (/c-i, fci) with —oo < 
/c_i < /cq < fci < 0. Moreover the size inf{\k^i — fcoM^o ^ ^il} of the time 
interval on which the solution is guaranteed to exist is controlled from below 
by l/i'{ko), ln|/co|, V{go,KQ) and £i(fco). 

2. There is continuity with respect to the initial conditions if we measure the 
space of initial conditions with the H\ x norm and the space of flow 
solutions with the BR-norm. 

3. Because of item 1 above, we have the following continuity principle: a flow 
solution ({g,K), {N,X)){k) is defined until past of k* < (or before fc, < i/ 
the flow is running in the past direction) iff limsupk^k-'^/'^ + ^{k) +£i{k) < 

00. 

Remark 10 i. As an outcome of the proof it is seen that if an initial state is C°° 
the flow solution is C°°. Moreover the growth of high-order Bel- Robinson energies 
Si (i > 2) is at most exponential (i.e. £i < Ce''*, z > 2) on any compact subinterval 
of I with constants C and r depending on the BR-norm of the solution on the 
subinterval. 

ii. The Theorem [5] does not make any a priori hypothesis on the topology of S. 
On manifolds with F(S) < we know the range of k is always a subset of (— oo, 0) 
but on manifolds with > the flow may be defined until ki — and beyond. 

The Theorem [5] does not discuss a continuity criteria in this case (at k = ki =0). 
One should be able to prove however that if ||fir||^2 > M > all along the flow and 
limfe_>fej=o l/^(fc) + V(^) + 'S'i(fc) < oo then the flow is extendible behind fci = 0. 
This situation corresponds to CMC solutions which are not time symmetric at the 
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slice {k — ki ^ 0} (i.e. K ^ Q at {k ^ 0}). The essential point in this case is that 
the lapse equation is solvable and the lapse N corresponding to the choice of time 
t = k is not degenerating under the assumptions. If the condition above on Hi^lj^,! 
is not imposed one may still be able to elaborate a criteria to guarantee that the 
flow is extended beyond fci = in which case the CMC solution is time-symmetric. 
In this case one should work with the choice of time t = k^ . Non of these issues 
we will treat in this article. 

iii. The Theorem[5]above should hold for initial states with W x regularity 
{i > 4) as well, namely they should give rise to a unique i?*-flow solution and the 
continuity criteria being the same as for x regularity. 

iv. A comment on the uniqueness of the solutions is in order. Call S the 
map that returns C°°-solutions from C°° initial states. Supply the C°° initial 
states with the Hj^ x iJji-norm and the C°° solutions with the BR-norm (for an 
appropriate time interval). As will be shown, the map S is Lipschitz and therefore 
admitting a unique extension to x H^. In this sense the solutions are uniquely 
determined from the initial data. But also, as will be proved (and as is required in 
the definition of flow solutions) , the flow solutions thus constructed are also space- 
time solutions. It follows from a well known result (see for instance |14j Theorem 
4.27 and references therein) that any space-time solution with the same initial data 
is diffeomorphic to the space-time solution arising from the flow solution. 



We will prove Theorem [5] towards the future, the proof towards the past is 
naturally the same. We choose t = k. 

Item 1. Step 1. Smooth the atlas A to make it C°° but compatible with A at 
the 4-th level of regularity. We prove that there are e, A and / such that for any 
C°° initial state {g,K) with II.9 — .9o||//^ < e and \\K ~ Ko\\ij2^ < e the C°°-solution 
is defined on / and has BR-norm bounded by A. For that we start proving 
that there are e, A and / such that while the C°°-flow solution is defined on /, the 
BR-norm over the interval where the solution is defined is bounded by A. Take A 

equal to C(||.go||H3^ + ||ifo||Hi + ||(£^o(io),So(to))||HO^ + 11(^^1 (io), Si (to))||//>;,) where 
C is a constant greater than one to be fixed later. Suppose there exists {gi, Ki) C°° 
initial states converging in x to (go, ^0) such that the BR-norm || (5, K)\\br 
over a time interval / = [to, ti] is equal to 2A, with ti t^. We can write 



Observe that the integrands are controlled by A by Propositions[29l Recall too that 
||J(Wi)||i2 and ||n|j//2^ are controlled by A. It follows from the Gauss equation 
that sup{te[to,ti]]{^i{t) — ^ 1(^0)} ^ as t^ ^ to- We see therefore that when ti — to 



Proof: 





to 
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is sufficiently small, for a constant C that was fixed big enouglo, the BR-norm of 
the solution on the interval [to, U] is less than A which is a contradiction. 

Next we prove that for such e, A and / the C°°-solution (/),g is defined in all of /. 
To show that it is enough to prove that Qi, i > 2 remains bounded (while defined) 
on /. In fact by Theorem[4]and Proposition[2]we know that: i. if the flow is defined 
until tt, in I, ii. its BR-norm is bounded by A and all Qi for i > remain controlled 
by A, then the flow {g,K){t) must converge in C°° to a C°° state. Therefore the 
flow can be continued beyond t* which is a contradiction. Using Proposition 1271 
bound the integrand involving J(Wi) in the Gauss equation for as 

JjNiEf3{W,Ub + Bf3iW,):^,)\dv, <Ci(£,_i,i.,||if||i.)Q?(Q? 

+ C2i£^-l,lyA\K\\L-)). 

Thus, from the Gauss equation we get the inequality 

\Q^■\<Cl{£,-l,A){Q, + C2{£^-l,A)Qh■ 

Proceeding inductively in i > 2 using this inequality we get that Qi for i > 2 has 
at most exponential growth (with constants controlled by A) as desired. 

Step 2: Until now we have a well defined map S from C°° initial states e-close 
to {go,Ko) in x Hj^ into C°° solutions defined on / with uniformly bounded 
BR-norm. We will prove now that the map S is Lipschitz if we take e sufficiently 
small. It will follow from the completeness of the space of flow solutions that S can 
be uniquely extended to a continuous map from x initial states e-close to 
{go,Ko) in Hj^ x Hj^ into H'^-How solutions in the sense of Definition [H To show 
that S is Lipschitz we proceed as follows. First we prove that if e is sufficiently 
small, any two states {gi,Ki) and (52, K2) in the ball of center {go, Kq) and radius 
e in H\ x H"^ can be joined by a path (,g, if)(A), with A G [0,1], of C°° initial 
states and with A-derivative controlled in H\ x H\ by || (50, -K'o)||H-]^x_ff2i ^^-^ 
two C°° states {gi,Ki) and (32,^2) as before, denote by g(A) the solutions with 
initial states (g, K){\), A € [0, 1]). Then we prove that we can choose e sufficiently 
small in such a way that there is a subinterval /' of / such that the BR-norm of 
g' = ^ on r (see Proposition [30] item 3) 

Il(0*g)'l|i3i?. = I|g'llc0(/',2)(ff^) + \\K'\\cO(rA)(H^) + ^\\{E'j,B'j)\\cO{rfi)(Hj,), 

is as small as we want. By Proposition[3D]this is enough to guarantee that the map 
S is LipschitS 

Let us start with the construction of the paths {g,K){X). Define {g,K){X) = 
X{gi, Ki) + {1 — X){g2, K2). Then, for every A find the transverse traceless part of K 

26c is chosen big enough to have £i{to < C{\\{Eo{to), Bo{to))\\[jO^ + || (£1 (to), Bi (io)) ll^^n^ )■ 
2^By Proposition 1301 when we differentiate the components of g(A) with respect to A we do not 

loose regularity. The claim that the map is Lipschitz follows essentially by integrating the fields 

along the path. 
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with respect to g and call it Ktt (see [10] )■ Now note that if Kq ^ then for e > 
sufficiently small it is Ktt (there is continuity in the York decomposition) and 
if = then for e sufficiently small g is conformally deformable to a metric of 
constant negative scalar curvature. Thus, according to |20j . the energy constraint 
is in all cases conformally solvable. It is direct to see that the path of initial 
states {g^K)[X) constructed in this way has || (g', iir')|j^3^x//^"iiorm controlled by 
the U\ X i/3i-norm of (go, ^^^o)- 

We next show that e can be chosen small enough in such a way that at every 
A G [0, 1] and for every path of C°° initial states as described above, the BR-norm 
of g' on a uniform interval /' C / is as small as we want. We will get control on 
llg'llsi? from a contradiction argument. We need some preliminary discussion on 
the A-derivative of the Weyl fields W — Wq and W = Wi . 

We first note that the A-derivative of a Weyl field is not necessarily a Weyl field 
as we have 

Note that {y^ ahcd)' = 0. We can construct from W a Weyl field by making 
it traceless according to the formula 

(121) W',,,, = Wl,,, - i(gacW^d + g^Wl, - gfc,Wl, - gadW^J, 

where we have defined WJ^^ = W^j^^^g'"', which makes it a (2, 0) symmetric and 
traceless tensor. We conclude from this that if we control g' in }i\ (i.e. all of 
its components at every time slice) and Q(W') then we control ||W'||^o^ (at every 
time slice). We will get control on Q(W') using the Gauss equation and to use the 
Gauss equations we need an expression for the divergence of W. We compute 

Jhcd = V"Wl,,, =V'^w;,,, - ^(gacV'^WL + gbdV^Wl, - gbcV"Wl, 
that we arrange in the form 

(122) J,,, = v^w;.,,, - i(v,w^, - v,w;j - i(g,,v''w;, - g^cV^w;,). 

We have therefore three different terms that we have to estimate on the RHS of 
the last equation. The first term can be computed by 

(123) v^w;,,,, = - v'-^w,,,, + J',,,. 

if we compute J^^^^. Contracting the last equation we get 

(124) v^w;., = j;,,,g^'^. 

which can be used instead of each summand in the third term of equation (|122p if 
we compute J^j^^. The second term VcWj^^^ — VdWj,^ in equation (|122|) is calculated 
by differentiating and contracting the identity 

1 *f 
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Differentiating we get 
and rearranging terms 



V„W;, - V.W;,, = - {V'^Wbede + V',W,ade + V,Wabde)S 

(125) , 1 f 1 f 

- y'^^Lea + ^^'fabc^*L + ^^fabeJ*'de 

Now, J (Wo)' = and J(Wi)' has an expression of the form 



(126) J(Wi)' = n'*VWo + n*V'Wo+n*VW^+r'*Rm*Wo + T*Rm*W^. 

Recah finaUy that J*^^ = ^JaZme'™c- Using the discussion above, let us give a 
schematic representation of the currents J(Wo) and J(W']^). The expressions will 
be used in the Gauss equations for Wq and W^. The expression for J(Wg) comes 
from using equations ()123p . ()124p and ()125p together with J(Wo)' = in equation 
(|122p . We can write (schematically) 

(127) J(W^,) = V'*Wo. 

To get a schematic expression for J(W'j^) use again equations (|123p . p24p and (|125p 
in equation (|122p . We can write (schematically) 

J(W;) = V' * Wi + J(Wi)' * g + e' * e * J(Wi) + e * e * J(Wi)'. 

with the expression for J(Wi)' borrowed from equation (|126p . 

From Proposition [30] (item 3) it is direct to show (see also iv. and v. in Remark 
[5]) that the BR-norms of g and g' over an interval /' control the norms 

||V'(t)||^^, ||n'(i)||Hi, l|W[,(t)||^^, ||e'(0||Hi, l|T'(t)|Uj. 

for every t in /'. Moreover, as Hg'llifjj is controlled by the BR-norms of g and g' 
(see V. in Remark [9]), we have 

(128) ||g'(t)-g'(to)||H^ <C(i-M, 

for every t in I' , where C depends on the BR- norm of g and g' on /'. The Gauss 
equation applied to Wq gives the inequality 

(129) |Q(W^)1 < CQ(W',), 

where C, again, depends on the BR-norms of g and g' on /'. To see that, use the 
schematic expression (jl27l) for the current J(Wq) in the Gauss equation (equation 

m) 
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together with the control on the norms ||n||^2^ and |1V'||^2^ mentioned before. 
Similarly one obtains an inequality for the time derivative of the Bel-Robinson 
energy associated to W'^ 

(130) IQ(W;)1 < C{Q{W[) + Q{W[)i). 

C, as above, depends on the BR-norms of g and g' over /'. Now, from the inequal- 
ities ([T29| and (fTSOl) we have 

(131) IQ(W^)(i) - Q(W^)(io)| < C(t - to), 

(132) IQ(W;)(i) - Q(W'i)(to)| < C{t - to). 

where, again, C depends on the BR-norms of g and g' on /'. Fix two A-independent 
horizontal vector fields V"' and . For any Weyl field W — Wq or W = Wi we 
have 

(133) Kfc^"^'' = ^'acbdT''T'^V''W^ + W acbdiT'^T'^VW'' + W acbdT''T"^V''W\ 
and recall (see iv. in Remark [O]) 

X"' N' 

(134) T'" = -— - — r». 
On the other hand 

(135) Wl,,, - Wabcd + ^(gacW^, + SmWI, - gfc,W;, - gadWiJ, 

and 

(136) W;, = -Wabcds""". 

From (fT33l) . (fm]) . (fTSS]! and (fTSe]) we can writeEl for E = E{Wq) or E ^ E{Wi) 
l|£^'ll?^o <C(A)(||W'||?,o + ||W||?,o||g'||?,0 

(137 3 3 A 

<C(A)(Q(W') + Q(W)||g'||^3^). 
Doing the same for B(Wo) or i?(Wi) and putting altogether gives 

\\{EiB'„)\\H<>^ + \\iE[,B[)\\HO^ <C{A){qHw',) + qHw[) + (Q^(Wo) 

(138) , 

+ Q^(Wi))||g'||«^). 

Note finally that by Proposition [30] |1 ^5' ||//2 and are controlled by the 

BR- norm of g and g' over an interval /'. This implies 

(139) \\9'{t)~g'{to)\\Hl <C{t-to), 



■^^Note that inside expression JTSTj there are norms with respect to A and norms with respect 
to g. We use implicitly that for any space-time tensor U we have yl(A)||U||£f2 < ||U||^2 < 
S(A)||U||„2. 
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(140) 



\\K'{t)^K'{to)\\Hl,<Cit~to), 



where C, as above, depends on the BR- norm of g and g' on /'. 

We start the argument by contradiction. We wiU end up showing that the BR- 
norm of g' over a time interval I' <Z I can be made as small as we want if e is 
chosen small enough. Suppose that there is i^i^i) and (g2,i, -^'2,1) converging 
to ((70) Kq) in Hj^ x i?^ and such that for some A.; G [0, 1] the BR- norm of g^(Ai) 
over an interval /' = [Iq, i,], with ti —^tg, is equal to a positive constant C. We note 
that as {gj_i,Kj^i), j = 1,2 converge to {go,Ko) in x it is {g[,K'^{\) 
in Hj^ X and uniformly in A € [0, 1]. This fact has several consequences on 
the A-derivative of the fields at time t — to. Namely, for the path {g,K){X) it is 
ll^'(*o)||/f3^ and ||X'(to)||ff^ -> 0. From (fm]) we conclude that ||T'(to)||ff3^ ^ 
and therefore ||g'(to)||H3j ~* ^' ^'^^^ decompose Wo(io) and Wi(io) into vertical 
and horizontal components using formulas (PT|) - (PS)) . Differentiating with respect 
to A we get that ||W[,(to)||^o^ -> and \\W[{to)\\H^^ ^ 0. It follows from (fmjl that 
W^oitoMlH'X and |!Wi(to)||HO 0. In particular we have Q(W^(to)) ^ 
and Q(W[{ta)) 0. Taking this into account, equations ^28]) . (fm]) . (fT32)) . (fT38)) 
and ()139p , (|140p , show that if U — to is taken small enough the BR-norm of g' over 
/' = [to, ti] is less than C which is a contradiction. 

Item 2. This item follows from the fact that the map S from initial states into 
flow solutions is Lipschitz. 

Item 3. This item follows from item 1. □ 



4.2 Long time flows. 

We give here an application showing that a (C°°) CMC- flow over a three- 
manifold with non-positive Yamabe invariant is a long-time flow if the correspond- 
ing globally hyperbolic space-time is future geodesically complete and the first order 
Bel-Robinson energy of the CMC flow remains bounded above. 

Theorem 6 Say F(E) < and (M ~ S x R, g) a smooth C°° maximally globally 
hyperbolic space-time. Say {g,K){k) with k e [a, 6) is a CMC flow where b is the 
lim sup of the range of k and say £1 < A. Then if (M, g) is future geodesically 
complete the CMC flow is a long time flow i.e. b ~ 0. 

A CMC slice with mean curvature ko will be denoted by Sfc^. Let t be a 
global smooth time function on M. Assume the range of t is (—00,00). Let T = 
— Vt/|g(Vt, Vi)| 2 be the unit future pointing normal to the foliation induced by 
t. Observe that if ^ : [c, d] -> M is a curve then i(^)' = dt{^) — g(Vt,^), therefore 
if g(T, ^) ^ everywhere the curve cannot be closed (i.e. ^(c) = S,{d)). 

Proof: We will proceed by contradiction and assume that 6 < 0. Under that 
assumption we first prove that the CMC foliation lies entirely between {t = to} 
and {t — ti}, for some to and ti. Suppose this is not the case, then there exists a 
sequence of points pi in E^. with ki ^ b and t(pi) 00. For a given i denote by 
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7i a past directed geodesic starting at pi and maximizing the Lorentzian distance 
to the initial CMC shce Sq. Clearly 7^ are perpendicular to Sq. Also because the 
lapse for the CMC time k is bounded above by 3/5^ the length of the geodesies 74 
is controlled above by h. There is a subsequence of geodesies ^i- converging into 
an inextendible geodesic 700 of bounded length and perpendicular to Eq. This is a 
contradiction as we have assumed the space-time was future geodesically complete. 

Now £1 is assumed less than A and if 6 < then FoZgi-j,-) (E) is bounded above 
for all k in [a, 6), therefore for the CMC flow not to be extendable beyond b there 
must be a sequence of points pi in E^. with ki ~* b such that the volume radius 
^g(ki){Pi) 0. Choose it in such a way that i^g{ki){p) > ^g(ki){Pi) for all p g 
Sfc.. We will be arguing with the rescaled space-time Afg with Ai = 
Denote with a subindex A^ any new rescaled quantity, for example k\. — ki/Xi, 
g\. = Xfg{ki) and so on. Noting that \K\'^dvg, \Ric\'^dvg and Qo scale as 
a distance~^ and Qi as a distance~^ , the sequence (T,,pi,g\^,K\^) is converging 
(strongly) in x into a flat (and complete) initial state (E,Poo, 5oo, -f^oo) 
where E is some three-manifold, goo is a flat metric and if 00 = 0. In particular 
there is sequence of geodesic loops U in E^^ (say of length one) and based at pt . 
Parameterized with the arc-length s we write U : [0,1] — > Eg^ C M. As noted 
above there must be for any i a, Si with gXi{li{si),TxJ — 0. At li{si) consider 
the orthonormal frame (^i, 62, 63, TaJ. Parallel transport it along li and call the 
resulting frames as es, 5 = 1,2,3,4. As observed above the states are becoming 
flattery therefore gA; (Vj .^i, e^) — > 0. Now think the loops k as inside the scaled 
space-time (M, Afg). Call ^oo(soo) a limit point of li{si) in M. As the scaled space- 
time is becoming closer and closer to Minkowski around Zoo(sco) the sequence of 
loops and frames converge to a solution of the system, i. gMinC^i loo,es) = 0, 
ii. V; es — 0, where gMin is the Minkowski metric on and loo '■ [0, 1] — > 
M*, and gMm(^(soo), K'Soo)) = 1- Necessarily, such solution must be a piece of 
segment of length one and therefore not closed. □ 

Remark 11 The same proof works to prove that a C°° CMC flow solution over a 
three-manifold with non-positive Yamabc invariant is cither a long time flow or the 
CMC foliation reaches the end of the maximally globally hyperbolic solution, in the 
sense that for any compact subset of the maximally globally hyperbolic space-time, 
there is a CMC slice which is disjoint from it (and to the future of a fixed CMC 
slice). 
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